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mergers 
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Abstract The merger of two neutron stars or of a neutron star and a black hole of- 
ten result in the ejection of a few percents of a solar mass of matter expanding at 
high speed in space. Being matter coming from the violent disruption of a neutron 
star, these ejecta are initially very dense, hot and extremely rich in neutrons. The 
few available protons form heavy nuclei (“seeds”) that absorb the more abundant 
free neutrons, increasing their size. The neutron density is so high that a substan- 
tial number of neutron captures occur before the resulting unstable nuclei can de- 
cay toward more stable configurations, converting neutrons into protons. Depending 
mostly on the initial neutron richness, this mechanism leads to the formation of up 
to half of the heavy elements that we observe in nature and it is called rapid neu- 
tron capture process (“r-process”). The prediction of the precise composition of the 
ejecta requires a detailed knowledge of the properties of very exotic nuclei, that 
have never been produced in a laboratory. Despite having long been a speculative 
scenario, nowadays several observational evidences point to compact binary merg- 
ers as one of the major sites where heavy elements are formed in the Universe. The 
most striking one was the detection of a kilonova following the merger of a neutron 
star binary: the light emitted by this astronomical transient is indeed powered by the 
radioactive decay of freshly synthesized neutron-rich nuclei and testifies the actual 
nature of compact binary mergers as cosmic forges. 
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Introduction 


After the discovery of the first stellar system formed by two neutron stars (NSs) 
[26], it was immediately realized that the decompression of NS matter following the 
merger of a binary neutron star (BNS) or of a black hole (BH)-neutron star (BHNS) 
system produces an ideal environment where the rapid neutron capture process (r- 
process) nucleosynthesis can take place. The very first calculations were carried out 
for BHNS binaries [34, 35]: in this scenario, during the last orbits of the gravita- 
tional wave (GW)-driven inspiral, the NS is tidally disrupted by the gravity of the 
more massive BH and a fraction of it is ejected into space. The coalescence of two 
NSs could also eject neutron-rich matter into the interstellar medium (ISM) through 
an even richer dynamics [12]. The first modeling of BNS merger nucleosynthesis 
was indeed accomplished a few years later [62, 17]. Since then, major progress has 
been made in understanding the mechanisms behind the ejection of matter and the 
properties of these ejecta. 

r-process nucleosynthesis is one of the fundamental processes responsible for the 
production of the heaviest elements in the Universe (e.g. [8, 9], see also [11] for a 
recent review). The binding energy per nucleon in nuclei increases almost steadily 
from lithium up to °°Fe and *°Ni. This allows the production of nuclei starting from 
H and He (mostly produced during the big bang) up to iron inside massive stars 
for increasing plasma temperatures in hydrostatic conditions (see e.g. [10]). The 
production of heavier elements through reactions involving charged nuclei would 
require even larger temperatures to overcome Coulomb repulsion. The kinetic en- 
ergy necessary to synthesise very heavy elements through fusion reactions becomes 
soon prohibitive for all plausible astrophysical scenarios. Moreover, at such high 
temperatures disintegration reactions would even dominate due to highly energetic 
photons. This strongly limits the nucleosynthesis yields produced through this path. 
The capture of a free neutron on a nucleus has instead the clear advantage of not hav- 
ing any Coulomb barrier to overcome and it is indeed the key process to produce the 
heaviest elements. However, free neutrons in non-degenerate conditions are unstable 
against B~ decay. Thus, neutron capture nucleosynthesis requires a source of neu- 
trons lasting for the relevant timescale over which the nucleosynthesis takes place. 
Moreover, the neutron-rich nuclei produced by neutron captures are unstable against 
B~ decay. While decay rates are constants, neutron capture rates crucially depend 
on the neutron density: if the neutron density is high enough (n, > 10°° cm~%), 
at least for a short timescale (~ 1s), neutrons are rapidly captured (increasing the 
mass number) before P decays increase the charge number and therefore produce 
isotopes of the next heavier element. Due to its rapidity, this process can happen 
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in explosive environments. This is the basic idea behind r-process nucleosynthesis 
and matter ejected from compact binary mergers provides precisely the conditions 
necessary for the r-process nucleosynthesis to occur. 

Over the past few years, several observational evidences have accumulated 
pointing to compact binary mergers as one of the main sites in the Universe 
where r-process nucleosynthesis takes place, including the observation of ele- 
ments synthesised through the r-process in the atmosphere of very old metal poor 
stars (see e.g. [61]). The first unambiguous detection of a kilonova (also called 
macronova), AT2017gfo, as one of the electromagnetic counterparts of the GW sig- 
nal GW 170817 (compatible with a BNS merger) represented the strongest evidence, 
so far, of this picture (e.g. [3, 60, 31]. 

In this chapter, we present the most relevant aspects of the r-process nucleosyn- 
thesis in compact binary mergers. We start by reviewing the conditions of matter 
expelled during a BNS or a BHNS merger. After that, we present how r-process nu- 
cleosynthesis proceeds in these ejecta. Finally, we overview the main observational 
evidences supporting compact binary mergers as major astrophysical sites for r- 
process nucleosynthesis. In most of our calculations we use cgs units. The physical 
constants employed through the text are the speed of light c, the reduced Planck 
constant fi, the Boltzmann constant kg, the gravitational constant G, the Stefan- 
Boltzmann constant Osp, the solar mass Mo, the masses of the electron, proton and 
neutron Me, Mp, Mn, and a generic baryon mass mp, which for our purposes can be 
assumed mp ~ mp. Temperatures are expressed both in Kelvin and in MeV (i.e. as 
kgT), depending on the context. We recall here that the conversion factor between 
GK and MeV is roughly one tenth, i.e. kg x 1GK ~ 0.086 MeV. The distribution of 
a generic nuclear species i can be expressed either in terms of its mass fraction X; or 
number abundance Y;. The former is defined as the ratio between the mass (density) 
of the species i over the total mass (density): X; = m;/mtot = P;/P, while the latter as 
the ratio between the number (density) of the i species over the total baryon number 
(density), Y; = N;/Np = ni/np. Clearly, np ~ p/m, and Y; = X;/A;, where A; is the 
atomic number of the species i. The electron abundance is defined as Y; = ne /np, 
where ne is the net density of electrons (i.e. the density of electrons minus the one of 
positrons). Due to charge neutrality, Ye = Yp free + Yp,nuclei: Where Y, free and Yp bound 
are the abundance of free protons and of protons bound in nuclei, respectively. 


Matter ejection from compact binary mergers 


Matter ejection in compact binary mergers happens through different channels. 
These channels are characterized by specific ejection mechanisms, which operate 
on different timescales and leave an imprint on the ejecta properties and, ultimately, 
on the nucleosynthesis. In the following, we review the properties of the ejecta from 
BNS and BHNS mergers by directly relating them to their merger dynamics. Before 
doing that, we briefly present some of the most relevant features that characterize the 
modeling of compact binary mergers and some of the more fundamental processes 
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that influence the ejecta properties. For more extended and complete information 
about these topics, we refer to the dedicated Chapters and to a few recent reviews, 
e.g. [59, 53], where detailed references to the original works can be found. 

The dynamics of the merger and of the ejecta expulsion depends on several in- 
trinsic parameters of the binary: first of all, on the nature of the coalescing objects 
(i.e. if it is a BNS or a BHNS system), but also on their masses and spins. Another 
relevant ingredient is the still uncertain nuclear equation of state (EOS) for matter at 
supranuclear density [46]. Any quantitative statement (and even a robust qualitative 
understanding) about the merger dynamics relies on detailed numerical simulations. 
The latter solve the equations of relativistic neutrino-radiation (v-radiation) hydro- 
dynamics coupled with dynamical space-time evolution. The hydrodynamics equa- 
tions are closed by a finite temperature, composition dependent EOS describing the 
microphysical properties of matter for a rest mass density that varies between stellar 
densities (a few g cm~*) up to several times 10!°g cm~?, corresponding to more 
than 10 times nuclear saturation density. 

Both in BNS and BHNS mergers an accretion disk around a central remnant is 
expected to form after the merger. In the case of BHNS systems, the remnant is al- 
ways represented by a BH, while in the case of BNS systems a massive NS, possibly 
collapsing to a BH on a variable timescale, usually forms. While the inspiraling NSs 
are in cold neutrino-less weak equilibrium, hot matter inside the merger remnant is 
out of equilibrium and its neutron-to-proton content changes due to neutrino-matter 
interactions. Several hydrodynamics processes increase matter temperature during 
the merger and the relevant temperatures range between 0 and ~150 MeV. Neutrino 
production is strongly boosted in hot and dense matter. Due to their low opacity, 
neutrinos become the dominant cooling source and their luminosity is of the order 
of a few 10% erg s~!, at least at merger and during the early aftermath. The decom- 
pression and heating of neutron-rich matter favors initially Ve luminosity through 
the reaction n+e* — p + Ve on the thermally produced positrons and Ly, > 2Ly,. 
The resulting net effect is to increase the electron fraction (leptonization), balanc- 
ing the proton-to-neutron ratio such that at later times Ly, ~ Ly,. In the densest part 
of the remnant, where matter density is above p ~ 10!?g cm~?, and the tempera- 
ture Trem is of the order of 10 MeV, the neutrino mean free path Zy = 1/(ngoy) 
is smaller than the size of the system (~ 10-100 km) for thermal neutrinos of en- 
ergy Ey ~ 3.15Trem : ly ~ 250 m(p/10!2g cm-3) | (Trem/10 MeV)~* . For this 
estimate we have used an approximated expression for the cross-section of neutrino 
scattering off free nucleons, oy © 69(Ey/m,c”)? with oo = 1.76 x 10-*4cm?, to 
evaluate the neutrino-matter cross-section. Absorption cross-section on free baryons 
have similar magnitudes and dependences. Then deep inside the remnant neutrinos 
equilibrate with matter and diffuse out on the diffusion timescale (~ seconds). Due 
to the presence of an accretion disk, neutrinos are emitted preferentially along the 
polar direction rather than along the equator, such that the polar flux can be a few 
times the equatorial one. The most relevant neutrino decoupling surfaces are located 
in the density interval 10!!~!*g cm~3, with v,’s decoupling at lower densities and 
temperatures than Ve’s and Vy ,r’s, since the neutron richness favors n+ Ve > p+e— 
as absorption process over other inelastic neutrino-matter processes. Typical neu- 
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Equal mass BNS, Mayn = 3.84 x 10° Mo 
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Fig. 1 Color coded two-dimensional histograms of the conditions of the dynamical ejecta, as ob- 
tained by an equal mass BNS merger simulation (BLh q = 1 run presented in [7]). On the left panel, 
the ejecta are characterized in terms of their specific entropy and electron fraction, while on the 
right panels of their specific entropy and expansion timescale. The top panels refer to an angular 
slice close to the rotational axis of the binary, while the bottom one to a slice close to the equator. 


trino mean energies in the decoupling region are Ey, ~ 10 MeV, Ey, ~ 15 MeV and 
Ey,,, ~ 20 MeV. Matter at even lower density, usually located at larger distances 
from the center, is irradiated by the neutrinos emitted at the inner decoupling sur- 
faces and this irradiation can change the neutron-to-proton content (i.e. Yo) through 
neutrino absorption on neutrons, protons, and nuclei. 

Finally, we recall that NSs are magnetized objects. During the merger and the 
subsequent remnant evolution, several mechanisms (e.g. dynamo amplification, 
magneto-rotational instabilities, Kelvin-Helmholtz instabilities) amplify the field 
strength. Even for initially low magnetic field (B ~ 10°~!°G), the field can be ampli- 
fied up to several ~ 10!5G and it becomes dynamically relevant during the merger 
aftermath, also for matter ejection. 


Ejecta from binary neutron stars mergers 


The first kind of ejecta emerging from a BNS merger are the dynamical ejecta. The 
dynamical timescale of the merger is set by the orbital (angular) velocity at the last 
orbit before merger, Vayn (Qayn), aS tayn ~ 27 /Qayn ~ (47% Rys) /Vayn Where Rys is 
the NS radius. Assuming a Keplerian behavior, vayn ~ \/GM/(2Rns), we obtain: 
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Unequal BNS, mayn = 6.81 x 104M 
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Fig. 2 Same as in Figure 1, but for a very unequal mass BNS merger (BLh q = 1.8 run presented 
in [7]). These conditions are qualitatively similar to the onse observed also in BHNS mergers. 
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for a typical total binary mass M of 2.7 Maand a NS radius of 12km [2, 4]. During 
the last orbits, when the two NSs approach each other, each of them gets deformed 
by the tidal field of the companion. As soon as the NSs touch, a large fraction of their 
kinetic orbital energy is converted into internal energy, while the tidal tails retain 
their orbital speed, Eq. (1). Since this is larger than the radial escape velocity, mat- 
ter in the tails is ballistically expelled with velocity vej + 0.1-0.3¢ S vayn. This is the 
tidal component of the dynamical ejecta and it mainly develops across the equato- 
rial plane of the binary. The temperature in these ejecta is only marginally increased 
by the tidal compression that precedes the merger (T < 1 MeV) and the emission 
of neutrinos is too weak to change Y, significantly. Thus, these ejecta mostly retain 
their original Y,, 0.05 < Y, S 0.15, and they increase their specific entropy only up 
to a few kg per baryon, s < 5kg baryon™!. If the total mass of the binary is too large 
for the matter pressure and for the rotational support to sustain the forming rem- 
nant, the latter collapses immediately to a BH. Otherwise, a (possibly metastable) 
massive NS forms in the center. This object is far from equilibrium and it bounces 
as the result of the gravitational pull and matter pressure response. As the sound 
waves generated by these oscillations travel through the remnant and reach lower 
density regions, they convert into shock waves, triggering the ejection of matter 


and 
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heated by compression and shocks from the outer edge of the remnant. These ejecta 
are called shock-heated ejecta and present a broad distribution of expansion veloc- 
ity, usually peaked around 0.2-0.3c, but with a possible high velocity tail extending 
up to 0.6-0.8c. Due to the action of shocks and compression the ejecta entropy in- 
creases, typically reaching values around 10-20 kg baryon”, but with a low-mass, 
high-entropy, high speed tail extending up to ~100 kg baryon ' The corresponding 
increase in temperature (initially, up to several tens of MeV before dropping due to 
matter expansion) produces a large density of electron-positron pairs and determines 
an increase in Y, in the ejecta due to positron captures on neutrons. Moreover, neu- 
trino irradiation coming from the forming remnant can further increase Y, through 
Ve absorption on neutrons [70]. Since neutrino emission is more efficient along the 
polar direction, the effect of irradiation is more evident at high latitudes. The com- 
bined effect of tidal tail interactions, hydrodynamics shocks and weak processes is 
the expansion of the dynamical ejecta over the entire solid angle, still with a prefer- 
ence along the equatorial plane (the mass distribution retains a sin? 0 dependence on 
the polar angle 0), with a clear gradient in the Y, distribution moving from the equa- 
tor (Y; ~ 0.1) to the poles (Y; < 0.4). In particular, matter above 0 ~ 45° is expected 
to have Y, = 0.25. The ejection of dynamical ejecta lasts for a few ms after merger, 
see Eq. (2), and its amount ranges between ~ 1074 and ~ 107? Mọ, depending on 
the binary properties and on the nuclear EOS (see e.g. [54]). The tidal component 
is more relevant if the high density part of the EOS is rather stiff or if the two NSs 
in the binary have very different masses. In this case, at least one of the two NSs is 
not very compact (i.e. Rys is larger) and the tidal disruption is very effective. The 
presence of shock-heated ejecta is instead more relevant in the case of equal mass 
mergers and for a soft nuclear EOS. In these cases, the NSs are more compact (i.e. 
Rys is smaller), the collision velocity at merger is larger (see Eq. 1) and the shocks 
are more violent. Typical conditions of the dynamical ejecta are presented in Fig- 
ures | and 2 for both an equal and a very unequal mass BNS merger, as obtained by 
detailed merger simulations in Numerical Relativity, at two different polar angles . 

While the massive NS forms in the center, matter compressed and heated up at 
the contact interface between the two NS cores is expelled outward. Conservation of 
angular momentum drives the formation of a rotationally supported, thick accretion 
disk of radial and vertical extension Rgisk and Haisk, such that its aspect ratio is 
(H/R)disk ~ 1/3. The typical disk mass Mgisk (where the disk is usually defined 
as the part of the remnant with density below 10!°g cm~? or, in the presence of a 
BH, outside of the horizon) ranges between 107° and 0.3M.. A prompt collapse 
to BH stops the disk formation, leading to lighter disks. An interesting exception is 
represented by very asymmetric binaries: in this case, the sudden BH formation is 
accompanied by a very efficient tidal disruption of the secondary NS, such that a 
significant fraction of it settles in Keplerian orbital motion outside the BH horizon 
[7]. 

The evolution of the rotating disk is governed by hydrodynamics, magnetic and 
weak processes. Several mechanisms are responsible for the local amplification of 
the magnetic field. The resulting viscosity of turbulent origin drives matter accretion 
onto the central object on the accretion timescale: 
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t oe H Ta x0 Ce H/R E Men a Raisk ae 
aeo aR a 0.02 1/3 2.6Mo 100 km 
(3) 


where œ is an effective viscosity parameter, Qgisk the Keplerian angular velocity, 
and for the remnant mass M,em and disk radial scale Rqisk we set as characteristic 
values 2.6Moand 100 km, respectively. 

During its secular evolution, the processes that determine the disk evolution 
can produce mass outflows, known as disk wind ejecta. Neutrino absorption, for 
example, redistributes energy and momentum inside the remnant, from hot and 
dense regions (p > 10!*g cm~*) to regions where the density decreases to p ~ 
10!0-!!g cm~’. This process inflates the disk, mainly in the vertical direction 
(neutrino-driven winds). At the same time, the accretion process implies an angu- 
lar momentum redistribution inside the remnant: while mass has a net inflow and 
the bulk of the disk is accreted, a fraction of it expands radially (viscosity-driven 
winds). In both cases, the disk expansion determines a drop in temperature and when 
T X 5-— 6GK free neutrons and protons recombine first in @ particles and then in 
heavier nuclei. The energy released is on average Enuc| © 8.6 MeV baryon”! and the 
corresponding expansion velocity can be estimated by equating the kinetic energy 
at infinity with the sum of the gravitational and nuclear energy released by recombi- 
nation: v/c S \/2 (Enuci/mpc2 — GMrem/(Rc?)), where Rrec is the radial scale where 
recombination occurs. Assuming Mrem = 2.6Mo, and Rrec © 450 — 600 km for disk 
winds, one obtains v ~ 0.034 — 0.074 c. Neutrino-driven winds can emerge on a 
few tens of millisecond timescale, while viscosity-driven winds on the longer accre- 
tion timescale. The earlier the wind develops, the faster the ejecta travel. Indeed, at 
earlier time the disk is hotter and the recombination radius is larger. 

Other processes happening inside the disk can also drive disk winds. For ex- 
ample, if the central remnant does not collapse to a BH, non-axisymmetric bars 
extending inside the disk (mainly m = 1 and m = 2 spiral modes) acts continuously 
on its innermost part, producing a net outflow of angular momentum that crosses 
the disk and expels matter from its edge (spiral-wave winds). Spiral-wave winds de- 
velops immediately after disk formation and possibly last up to the point where the 
central massive NS collapses or the bars are dissipated by GW emission. Moreover, 
if in addition to disordered local fields, large scale magnetic fields develop inside 
the remnant, magnetic pressure and the Lorentz force can further accelerate matter 
producing magnetically-driven wind disks. Both in the magnetic and in the spiral- 
wave wind cases, expansion velocities are intermediate between the fast dynamical 
ejecta and the slower recombination disk winds, v ~ 0.1 —0.2c. 

The many processes taking place inside the disk are very effective in unbinding 
mass from it. The scale that sets the ejecta amount is the mass of the disk itself. 
While neutrinos alone are able to unbind only a few percents of the disk (and even 
less in absence of the very luminous massive NS), the other mechanisms unbind 
between 0.1 and 0.4 Maisk. It is worth mentioning that these mechanisms can work 
at the same time, with disk consumption being the only really competing factor. 
Disk-wind ejecta are thus likely the most relevant source of ejecta for BNS mergers, 
as visible in Figure 3. 
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Fig. 3 Comparison between the dynamical and the disk-wind ejecta, as obtained by a large set of 
BNS merger models, employing several nuclear EOSs. While the amount of dynamical ejecta are 
computed within the simulations, the mass of the wind-disk is assumed to be 20% of the mass of 
the accretion disk at the end of the simulations. Figure taken from [54]. 


The ejection timescale of disk winds is comparable to the weak reaction timescale. 
As for the dynamical ejecta, while the initial Y, in the disk is set by the cold weak 
equilibrium of the merging NSs, the hot temperature increases Y, inside the ex- 
panding winds due to positron and neutrino absorption on neutrons. Assuming that 
neutrino irradiation is effective and long enough to reach equilibrium, Y, tends pro- 
gressively toward [52, 42]: 


Ly, Wy, Ev, — 2A +1.2A2/ey, ) a r 


Yera © (: + Ly, Wy, Ev, +2A + 1.2A2/éy, 
where A = (m, — mp)? ~ 1.29 MeV, & is the ratio between the average squared 
neutrino energy and the average neutrino energy (which for relevant spectral dis- 
tribution gives £y ~ 1.2 Ey), and W; ~ 1+7(E;)/mpc? with ny, = 1.01 and 
Ny, = —7.22 is the weak magnetism correction factor. For typical neutrino lumi- 
nosities and mean energies, Ye eq © 0.45 < 0.5. Thus, depending on the ejection 
time and on the strength of the neutrino irradiation, the electron fraction in the wind 
ejecta shows a broad distribution between 0.1 and 0.45. In the case of short-lived 
massive NS the bulk of the disk wind ejecta have Y, between 0.2 and 0.3, while 
in the case of long-lived remnant, the higher neutrino luminosity drives the Y, dis- 
tribution toward Y, = 0.3. Matter in the disk is shocked by waves produced by the 


~N 


central remnant and it is heated by viscous dissation. Then its entropy increases to 
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an average value of 15-25 kg baryon”, with possibly high entropy tails extending 


also in this case up to ~100 kg baryon™!. 


Ejecta from neutron star-black hole mergers 


The ejection of matter from BHNS systems shares many similarities with the one 
from BNS mergers, but it also present crucial differences. Astrophysical BHs are 
characterized by their mass, Mgy, and spin, Jgy. In the following we assume that 
the orbital angular momentum of the binary points toward the positive z direction 
and that the BH is significantly more massive than the NS (Mpu Z 5 Mo, while 
Mys <2 Mo). Assuming Jpy to be also along the z axis, we characterize it through 
the dimensionless spin parameter agy = +Jpuc/ (GM), where 0 <a < 1 if the 
spin points toward the positive z direction, —1 < a < 0 otherwise. When at the 
end of the GW-driven inspiral phase the two compact objects approach each other, 
the merger fate depends on the location of the BH last stable circular orbit, Rjsco 
(i.e. the radius inside of which circular time-like test-mass orbits in the equatorial 
plane become unstable to small perturbations), with respect to the tidal distance, 
didal (i.e. the BH-NS distance at which the gravitational force on a test mass at the 
NS surface equals the tidal force pulling the mass toward the BH). For Kerr BH 
Risco = Mpu f(a), where f(a) is a monotonically decreasing function such that 
1 < f(a) <9 and f(0) = 6, while didai ~ (2Mpu/Myns)! Rys. If didai < Risco 
the NS is swallowed by the BH before a significant tidal disruption can occur. No 
mass is practically left outside the horizon to form a disk or to become unbound. 
Otherwise the NS experiences a partial tidal disruption before most of its mass gets 
inside the BH horizon. Mass coming from the farther NS edge can become unbound 
in the form of tidal dynamical ejecta inside a crescent. The rest of the unswallowed 
mass sets into a Keplerian accretion torus around a spinning BH. From the above 
relations, it is clear that the probability of leaving mass outside the BH horizon 
increases for stiffer nuclear EOS, smaller mass ratios, and larger agy . Numerical 
simulations show that 0 < Mayn < 0.15Mofor the dynamical ejecta while the mass in 
the torus is such that 0 < Mgisk S 0.5Mo. While the dynamical ejecta mass decreases 
for larger NS masses, the most massive torii are observed for very massive NSs 
(and for smaller mass ratios and larger spins). The NS gets compressed during the 
inspiral. However, the tidal nature of the merger keeps the entropy in the dynamical 
ejecta low (s ~ a few kg) such that kgT < 1 MeV always. The electron fraction of the 
ejecta stays also very close to its cold neutrino-less weak equilibrium value, i.e. Ye ~ 
0.05. Inside the torus, even in absence of shocks produced by the central remnant, 
accretion and disk dynamics heat up matter up to kgT ~ 10 MeV, producing ~ 
10°*erg s~! of accretion-powered neutrino luminosities. The subsequent production 
of disk wind ejecta is similar to the one observed in BNS mergers. Due to the lower 
neutrino luminosity, the electron fraction of the disk wind ejecta shows a broad 
distribution, 0.1 < Y; < 0.35, with a lower average value, compared with the wind 
disk ejecta produced by long-lived NS remnants. 
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Ejecta expansion and thermodynamics 


As we have seen both for BNS and for BHNS systems, during the merger cold 
and neutron-rich nuclear matter is heated up and leptonized by several processes, 
including matter compression, hydrodynamics shocks and neutrino irradiation. The 
bulk of the dynamical ejecta and of the remnant disk originate from inside the outer 
NS core, where p > 10!4g cm~? and Y, < 0.1. In these conditions, nuclei are fully 
dissociated in free neutrons and protons (homogeneous nuclear matter), and due to 
charge neutrality their initial abundances are related to the electron fraction by Y, = 
Y. S 0.1 and Y, = (1 — Ye) 2 0.9. After having reached its peak temperature, matter 
is expelled and the density and temperature drop, while neutrons and protons start 
to form nuclei. In the following we will always assume that the peak temperature of 
the ejecta exceeds 4 GK. 

Matter and radiation in the expanding ejecta can be considered as a neutrally 
charged plasma consisting of nuclei (often distinguished in free neutrons n’s, free 
protons p’s, Œœ particles and a distribution of all other nuclei), electrons, photons and 
neutrinos. During most of the relevant nucleosynthesis timescale, photons, electrons 
and nuclei are in thermodynamics equilibrium. It is thus useful to consider an ex- 
panding fluid element as a Lagrangian particle characterized by an evolving density 
p = p(t) and temperature T = T (t). Nuclei and nucleons can be treated as an ideal, 
non relativistic classical gas, obeying Maxwell-Boltzmann statistics. Electrons are 
degenerate in the early phase of the expansion before entering a classical gas phase, 
and they can be described by an ideal Fermi gas of arbitrary degeneracy. Photons 
are always characterized by a black body spectrum, typical of massless bosons in 
thermal equilibrium with matter. Once the density, temperature and composition are 
given, the entropy of the system can also be computed through the resulting equation 
of state, s = s(t). 

In the first phase of the ejecta expansion, when hydrodynamics processes are 
still active, the fluid density decreases approximately in an exponential way: p(t) = 
Po exp(—(t—to)/T), where po is the density at the onset of the expansion, t = fo. 
The expansion timescale, T, quantifies how fast the density drops during the first 
phase. Faster dynamical ejecta have smaller t’s (1 < t[ms] < 10) than slower disk 
winds (10 < t[ms] < 100). The expansion timescale can be related to the expansion 
velocity vej through: 


R R Vej \ T! 
~ — ~ 16.7 ms| o ) (È) 
sa S ms (n) 01c ©) 


where R is the lengthscale where matter becomes unbound, depending on the kind 
of ejecta and ejection mechanism. For the dynamical ejecta, the more impulsive 
expulsion set a lower R ~ 300 km than for the disk wind ejecta, R ~ 600 km. After 
a timescale thom ~ 3T, internal dynamics and momentum redistribution cease and 
fluid elements start to expand with an approximately constant velocity: the fastest 
fluid elements at the front, the slowest at the bottom. Such expansion profile is said 
homologous since v œ r and in this phase the density evolves as: 
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P (£) = Phom Ea SDE er (, x en | (az) (7) ` (6) 


with Phom Varying inside the range 104 — 108g cm7?, and the larger values are for 
low entropy, fast expanding ejecta. This density evolution is a very relevant input 
for the r-process [58]. After the nucleosynthesis, the internal energy of the ejecta 
is subdominant with respect to the kinetic energy and the expansion proceeds adi- 
abatically. Moreover, for the relevant density and temperature conditions, the EOS 
is dominated by the relativistic photons and electrons. For such a gas, an adiabatic 
expansion satisfies pT? ~ const and thus T(t) « t~!. While this temperature evo- 
lution profile is not accurate during the nucleosynthesis epoch, it provides a good 
approximation after it. 


r-process nucleosynthesis in compact mergers 


We now move to the study of the nucleosynthesis happening inside the ejecta. Be- 
fore speaking about the r-process, we introduce a few basics, although necessary, 
concepts in nuclear reaction theory. We address the reader to Refs. [10, 27, 41] for 
general introductions and to [67, 11] for more specific reviews on the r-process. 

Nuclear abundances inside the ejecta evolve in time as a result of nuclear (both 
strong and weak) and electromagnetic reactions. The usage of abundances in nu- 
cleosynthesis calculations is very sensible because, being the ratio of densities, it 
allows to decouple the effects of reactions on the composition from the effects due 
to matter expansion. Strong and electromagnetic reactions conserve separately the 
number of protons and neutrons. Strong nuclear reactions include fusion reactions 
among nuclei. In the case of transfer reactions, they are often indicated as B(i,o)C 
to emphasise the transfer particles i and o. For example, in a (n, œ) reaction a free 
neutron is absorbed by a nucleus B = (A, Z), and a nucleus C = (A—3,Z—2) and an 
œ particle are produced in the final state. Strong nuclear interactions include also œ 
decays, (A,Z) — (A — 4,Z — 2) + æ, and spontaneous or induced fission processes. 
The most relevant electromagnetic reactions involving nuclei are the photodisin- 
tegration reactions, B(y,o)C, and their inverse absorption processes, B(i, y)C. For 
example, in the case of a (y,n) reaction, the absorption of a photon on a nucleus 
B = (A,Z) produces a nucleus C = (A — 1, Z) together with a free neutron. The op- 
posite reaction (n,7), called neutron capture, consists in the absorption of a free 
neutron on a nucleus B = (A, Z), producing a nucleus C = (A+ 1,Z) and a photon. 

Weak nuclear reactions include B= decays, electron and neutrino captures, and 
they convert neutrons into protons and vice-versa, changing the electron fraction 
of matter. In the expanding neutron-rich ejecta, B~ decays, (A,Z) > (A,Z+1)+ 
e~ + Ve, are the most important weak reactions while in the initial expansion phase 
neutrino irradiation acts through neutrino absorption, as previously discussed. 

Any of the transfer and absorption reactions B(i,o)C is characterized by its cross- 
section, Og(;,5)c, Which in general depends on the energy of the colliding particles. In 
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an astrophysical plasma, the distribution of the colliding energy depends on the local 
thermodynamics properties and the reaction rate, rg(;,.)c (defined as the number of 
reactions occurring per unit time, per units volume and per reactant pair), is given 
by rg(i.0),c = (OY) Biio)c NBN: Where (OV) g(;,0)c is the product of the reaction cross 
section times the relative velocity between B and i, averaged over their thermal 
distributions. 

In the case of photon or neutrino absorptions on a nucleus B, the reaction rate 
can be expressed in terms of an effective destruction/decay rate, AB(y/v,0)C> as 
'B(y/v,0)C = AB(y/v.o)c MB, Where Agy is a function of T only (since the photon 
gas properties depend only on temperature), while Ag. depends on the neutrino 
spectrum, which is in general not in equilibrium with the plasma. Then, strong and 
electromagnetic reactions depend only on the local plasma properties. Alpha and 
beta decays are instead characterized by a constant decay rate Ag /p> related to the 
nucleus half-life t /2 by ha p= 2/ti /2- The lifetime of a nucleus against a certain 
reaction can be defined as the inverse of the corresponding rate, T; = 1/Aj. 

For large temperatures, strong and photodisintegration reactions are character- 
ized by large reaction rates. For T > 4—5 GK the resulting fusion and photodisin- 
tegration timescales become much shorter than the weak and dynamical timescales 
in the ejecta. Thus, these reactions can be considered in equilibrium among them 
and with their inverse reactions. This condition is called Nuclear Statistical Equi- 
librium (NSE) and the nuclear abundance in NSE are fully determined by the local 
thermodynamical conditions, i.e. by p, T and Y.. While p and T vary due to the 
expansion, Ye changes due to weak processes, but in both cases on much longer 
timescales than the nuclear NSE timescale. When the temperature decreases below 
~ 4—5 GK some reactions characterized by small Q-values become slow enough 
that NSE is no more guaranteed across the entire nuclear distribution, and especially 
for nuclei characterized by magic nuclear numbers, i.e. close to shell closure con- 
ditions. This transition is called NSE freeze-out. Since a large fraction of direct and 
inverse nuclear reactions are still very fast, the nuclear distribution splits into areas 
of Quasi Statistical Equilibrium (QSE), where equilibrium conditions still apply on 
sub-sets of nuclei. When the temperature decreases even further, no equilibrium ar- 
guments apply and fully out-of-equilibrium nucleosynthesis occurs. While in NSE 
conditions accurate abundances can be computed even without precise information 
on the reaction rates, a detailed knowledge of the properties of the ejecta and of 
nuclei all over the nuclear chart are requested to predict accurate abundances after 
NSE freeze-out. Nevertheless, equilibrium arguments still provide a useful tool to 
understand the basic feature of nucleosynthesis in QSE conditions. 


Compact binary mergers as r-process site 


Before entering the details of the r-process, we first motivate why the ejecta of 
compact binary mergers represent a suitable environment for r-process nucleosyn- 
theshis. The ejecta come from high density conditions and usually experience high 
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enough temperatures such that matter is mostly dissociated into free neutrons and 
protons under NSE conditions at their peak temperature. As temperature drops, neu- 
trons and protons recombine first into œ particles. The subsequent building of the 
most tightly bound iron group nuclei depends on the three body reactions responsi- 
ble for the assembly of heavier nuclei, namely 2a +n — °Be+yand 3a C+ y. 
The first one is the most relevant in neutron-rich conditions. Triple reactions are in 
competition with their inverse photodestruction reactions. Due to their triple nature, 
the former are favored by larger densities, while the latter are strongly enhanced by 
higher temperatures since for the photon density and mean energy one has ny œ p? 
and (Ey) « T, respectively. If the plasma is radiation-dominated, the value of the 
specific entropy ultimately determine whether iron group nuclei can form. Indeed, 
if s Z sy+s,+ © 2sy, the density can be computed as: 


4n?mp kT? y 6 3 T sS md 
ps2 (am) ie = 3.02 x 10° g cm (<x) (oz) ; (7) 
In very neutron rich conditions (Ye ~ 0.1), when temperature decreases between 5 
and 2.5 GK, œæn reactions occur more efficiently than their inverse photodestruc- 
tion reactions only for p > 3 x 10° g cm~?, i.e. for s < 100kg baryon~!. This implies 
that in low and moderate entropy conditions (typical of the merger ejecta) iron group 
nuclei are formed in NSE conditions, while in the high entropy tail œ particles are 
mostly produced when NSE equilibrium is no more guaranteed (@-rich freeze out). 
Almost all protons are bound inside nuclei and a distribution of heavy nuclei 
around the iron group (called seed nuclei) is present at the NSE freeze-out. For the 
seed nuclei, one can introduce a representative average nucleus ((A) seed, (Z) seed) 
defined such that: 


(A) seed = ( X. A faz) /Ysced ; (Z) seed = ( y Z Yaz) /Yseed ; (8) 
(A >4 


Z), A>4 (A,Z),A 


where we have defined Yseeq = Yiia.z) a>4 Yaz). While for Ye ~ 0.5 the nuclear 


abundance distribution has its peak around *°Fe and *°Ni (the most bound, sym- 
metric nuclei, characterized by the presence of a Z = 28 proton shell closure), for 
Y, < 0.5 the distribution moves toward more exotic, neutron-rich iron group nuclei. 
However, there is a limit for the amount of neutrons that can be bound inside a nu- 
cleus: for 20 < Z < 40, nuclei located at the neutron drip line have (Z/A) min ~ 0.3. 
Moreover, due to the neutron richness it becomes energetically favorable to have a 
certain fraction of free neutrons even for Y, > 0.3. Starting from baryon number con- 
servation written as X; + Xn seed + Xp, seeda © 1 (where X, /p,seed indicate the fractions 
of neutrons or protons bound in seed nuclei), from the charge neutrality condition 
Ye ~ L(A Z) A>4 (Z Y az) and using the above definitions, the free neutron fraction 
can be evaluated as: 


Yn at= ((A) seea/(Z) seed) Ye - (9) 
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Assuming for simplicity that (Z) seq 28, and that the dependence of (A) scea on Ye 
is approximately linear between A/Z = 0.5 and A/Z = (A/Z) min: 

(A) seca = 56 +70 (0.5 — Yo) , (10) 
we obtain an approximation for Y;,: 


Y, = 1=3.25Y,4- 2.52, (11) 


that varies between 0 for Y, = 0.5 and 1 for Y, = 0. Since almost all protons are 
contained inside the seed nuclei, we can also estimate the seed abundance as: 


Yeeed © Yo/(Z) seed © Ye /28. (12) 


Neutron capture is the key reaction to produce heavy nuclei beyond the iron 
group, once NSE freeze-out has occurred. In accordance with the definition of reac- 
tion rate and in analogy to the effective rate definition, the lifetime of a generic seed 
nucleus against a neutron capture reaction can be estimated as: 


Tny) = 1/ (ma(Ov) (ny) > (13) 


where n, is the free neutron density. The Q-value of an (n, y) reaction is the energy 
gained by the nucleus (A,Z) by acquiring a neutron and it is computed as Q = 
(maz + myc? —ma+41,z)) c?. It is also equal to the energy required to remove a 
neutron from a (A+ 1,Z) nucleus, called neutron separation energy of the (A+ 1,Z) 
nucleus, S;,(4+1,z). Neutron captures have Q-values that range from 0 for exotic 
neutron-rich nuclei at the neutron drip line up to ~ 15 MeV close to the valley 
of stability. The leading contribution to 0(,y) is provided by the s-wave term of 
the partial wave expansion: On ~ (x 7. k?) Ts, where T; ~ 4k/k' is the transmission 
coefficient obtained by considering a neutron moving against the potential barrier 
of the nucleus, and k,k’ are the wave numbers of the particles in the initial and final 
state, respectively. If u ~ m, is the reduced mass of the parent state and E ~ kgT the 
thermal energy scale, for non-relativistic energies (kgT ~ 0.1 MeV < Q < mc?) 
the relative speed is v = ,/2E/, and the wave numbers are k ~ \/2uUE/h and 


k ~ \/2u(S,+E)/h ~ \/2US,,/h, so that the typical cross-section is: 


keT yey g a 
~ 3.70 barn 2 14 
Sng) ~ 3-70 am (ey) 5 MeV a) 


Since 0, « 1//T and v = VT, (OV) ny is approximately constant for thermal, non- 
relativistic neutrons. Thus, the lifetime against neutron capture, Eq. (13), depends 
mainly on n,, which for a plasma of density p and free neutron fraction Y„ is simply 
nn = PYn/mp & 2.99 x 10°”^cm™? (p/10 g cm~’) (¥,/0.5), and T(n,y) becomes: 


= -1 1/2 
p Yn Sn 
wü ' 1 
Fan ~o s(a) (5) (z) ea 
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The density in the ejecta changes considerably with time as a result of the homol- 


ogous expansion, Eq (6). If we further consider that Phom ~ 10°g cm~? at NSE 
freeze-out for s ~ 10 kg baryon~!, see Eq (7), then the free neutron density at time 


t is: 
Y, T \37t \-3 
~ 8.08 x 10% cm~? ( —_Phom 16 
k j oe: (ta 0.5 Cre (; 2) i (16) 


and Y 


y) can be also expressed as a function of time after merger as: 


n, 


Yn Phom L T = Sn oe Bo 
~ 11 17 
Tiny) 8 us (= pen — | (a =) 5 MeV Gj ;) a7) 


The capture of one or more neutrons increases the mass number by one or more 
units without increasing the atomic number. Then, neutron captures move the nu- 
clear abundances toward the neutron-rich side of the nuclear chart. Two kinds of 
reactions compete with neutron capture in producing heavier and heavier neutron- 
rich nuclei: (7,7) reactions and the B~ decays. 

For the (7,7) reactions, while the high neutron density guarantees high (n,7) 
rates, high energy photons are required to knock a neutron off a nucleus, over- 
coming the neutron separation energy S,, and boosting the photodestruction rates. 
For T = 4GK the two set of reactions are in NSE, meaning that the temperature 
is large enough to bring (n,y)-(n,y) at equilibrium everywhere, also close to the 
valley of stability where S, ~ 8 — 10 MeV. In neutron-rich conditions, when the 
nuclear distribution is shifted toward the neutron-drip line, the relevant S,, can be as 
low as ~ 1 — 2 MeV. Thus, temperatures lower than 4 GK are enough to preserve 
the (n, y)-(n, y) equilibrium after NSE-freeze-out on the neutron rich side of the nu- 
clear chart. This is indeed the QSE typical of neutron-rich conditions. Since (y,n) 
reactions are the inverse of (n, y) reactions, their rates are related by detailed bal- 
ance conditions: (4,7) (OV) (n,y) =(A+1,z)A(y,n)> Where À(y n) is the photodisinte- 
gration rate of (y,n) reactions. Assuming these reactions to be in equilibrium, i.e. 
(A,Z) +n © (A+1,Z)+y, the chemical potential of the different nuclear species 
involved are related by U4 z) + Un = M(a+1,z) (we recall that Uy = 0). Using the 
expression of the chemical potentials for an ideal Maxwell-Boltzmann gas in its 
relativistic version (i.e., including the rest mass contribution), one obtains: 


nayz) — ( 2ah 3/2 Ga+ız(T) (A+1 ee Sn(A+1,Z) (18) 
N(A,Z) Nn E mpkgT 2 Gaz (T) A P kpT f 


where Gr, 7) (T) is the nuclear partition function, dependent on the matter tempera- 
ture. Assuming G(4+1,z)/Gya.z) ~ 1 and (A+ 1)/A ~ 1, we finally obtain an expres- 
sion for Tn) = 1 / Ayn)» the lifetime of a seed nucleus against photodestruction: 


-3/2 1/2 
kgT ) ( Sn ) exp (Sp/kpT ) (19) 


~ 0.19 us ( —<B- _ 
Tyn) j s (suey 5MeV exp(30) 
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where the reference temperature (corresponding to T ~ 2 GK) and neutron separa- 
tion energy are chosen such that S„/kgT ~ 30. This timescale depends heavily on 
the temperature due to the presence of the exponential factor that rises by orders of 
magnitudes as the temperature decreases. Additionally, it depends also on the nu- 
clear masses through S, and the latter sets the scale in the exponential argument, 
meaning that for a given temperature Ty) changes dramatically between the valley 
of stability and the drip line. 

Initially, at NSE freeze out (T ~ 4 GK), the temperature and the density are large 
enough such that both Tn,y and Ty n are much smaller than the dynamical timescale. 
The exponential term in Eq. (19) set the typical S, necessary to guarantee the equi- 
librium between the two reaction sets. As T and p decrease, Tyn) changes much 
more sensibly and (n, y)-(y,n) equilibrium can establish only down to T ~ 1 GK, 
assuming typical S, ~ 1-3 MeV. 

For the B~ decays, using the low-energy limit of the weak interaction theory at 
leading order, the transition matrix element .@ is proportional to the Fermi coupling 
constant, Gr. Remembering that in natural units Gp is the reciprocal of an energy 
squared, that the time is the reciprocal of an energy, and that the only relevant energy 
scale in the process is the Q-value of the decay, the lifetime of a nucleus against p 
decay, Tg, must be proportional to Q->. This dependence can be seen also as a 
consequence of the three-body nature of the final state. The Q-value of B~ decays 
involving neutron-rich nuclei is roughly proportional to the neutron excess, D = 
N —Z, and varies between a fraction of MeV close to the valley of stability and 
~5-15 MeV at the neutron drip line (with larger values at lower mass numbers, 
where the neutron excess can also be much larger). Using the decay of the free n as 
representative B~ reaction (for which Q = A and D = 1), we estimate the typical 
B-decay lifetime as: 


A 5 Q -5 D =5 


where in the last step we have further assumed that Q ~ DA. Then, B~ decays act 
on much longer timescale than neutron captures at NSE freeze-out and during the 
(n, y) — (y,n) equilibrium, and become competitive only when the temperature and 
density have significantly dropped. This difference in the neutron capture and in the 
B~ decay timescale qualifies the the ejecta from compact binary mergers as one of 
the astrophysical sites for r-process nucleosynthesis in the Universe. 

Starting from a certain seed nuclei distribution, the fast neutron and photon cap- 
tures move abundances within the same isotopic chain (i.e., the sequence of nuclei 
characterized by the same Z and by an increasing N). However, nuclear decay is a 
stochastic process and a fraction of nuclei decay from one isotopic chain to the next 
one even if Tg is significantly larger than 1, y), especially if the (n, y)-(y,7) equi- 
librium is maintained on a timescale comparable to or longer than tg. Moreover, Sn 
is not a smooth, monotonic function inside the nuclear chart. The closure of neutron 
shells inside the nucleus at the magic numbers N = 28, 50,82, 126 determines an in- 
crease of S, around those values. The corresponding nuclei become waiting points, 
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where T(n y) increases and matter tends to accumulate. For these nuclei the B~ decay 
starts earlier to be competitive and matter flows through them from an isotopic chain 
to the next one. 

The net results of the combination of neutron captures, photodestructions and 
p~ -decays on the whole distribution of nuclei emerging from NSE during the (n, /)- 
(y,n) equilibrium is a characteristic nuclear distribution that proceeds as a river 
inside the neutron-rich side of the nuclear chart. This is called the r-process path. 

The final point of the r-process nucleosynthesis depends on how many free neu- 
trons are available to be captured by the seed nuclei. This number is called the 
neutron-to-seed ratio, Y,/Yseeq, and the end point of the r-process can be estimated 
as: 

(A) final i (A) seed al Yn /Yseed , (21) 


where (AYseeq is the average mass number of the seed nuclei. For example, assum- 
ing (A)seea ~ 80 — 90 (see Eq. 10), to produce element with A ~ 130 it is necessary 
to have Y;,/Yseeq © 40 — 50; for A ~ 195, Yn/Yseea ~ 95 — 115 (we will see later 
than these are the mass numbers of the so-called second and third r-process peaks); 
while for uranium and thorium (A ~ 235) Yn /Yseea ~ 145 — 155. For low or moder- 
ate entropy ejecta, the value of Y, primary determines the neutron-to-seed ratio and, 
from that, how far the r-process nucleosynthesis proceeds in producing heavy ele- 
ments starting from iron group seed nuclei. Using Eqs. (11) and (12) we can obtain 
a simple estimate for the neutron-to-seed ratio as a function of Y, in low entropy 
conditions: 


Yn/Yseed ~ (Z)sced/Ye — (A) seed ~ 28/Ye — 70 Ye — 21. (22) 


This expression diverges for Ye —> 0, as there are no seeds, while it goes to 0 
for Y, = 0.5. For intermediate values Y, = 0.1,0.2,0.3,0.4 we obtain Y, /Yseed © 
250,105,51,21, respectively. These formulae, relying on the simple assumption 
of a linear dependence of (A)sceq on Ye, Eq (10), must be understood as very 
rough estimates that nevertheless catch the most relevant trends for typical entropy 
(s~ 10 kgbaryon~!) and expansion timescale (t ~ 10 ms) in the ejecta. More de- 
tailed and physically motivated calculations (e.g. [21, 38]) extending to broader 
ranges of possible conditions, show that the neutron-to-seed ratio has a more com- 
plex dependence on the (thermo)dynamics conditions at NSE freeze-out, in particu- 
lar, Yp /Yseea is larger for lower Y., smaller T and larger s. For example, for Ye = 0.4 
the neutron-to-seed ratio is such that ¥;,/Yoeea ~= 5° /¥3T. 


The working of the r-process in compact binary mergers 


The evolution of the nuclear abundances and the calculation of the nuclear energy 
released during the nucleosynthesis is computed by nuclear reaction networks. A 
nuclear network is a large system of coupled ordinary differential equations. For 
each nucleus i = (A;,Z;), its abundance Y; = Y(A;, Z;) evolves according to: 
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Fig. 4 Evolution of a few selected abundances computed for a fluid elements expanding in space 
with homologous expansion, see Eq. (6), and initially characterized by a specific entropy s = 
10 kgbaryon™!, an initial Y, = 0.10 and expansion timescale T = 10ms, corresponding to typical 
compact binary merger ejecta conditions (left panel). Heavy nuclei are defined as nuclei for which 
A > 120. In addition, relevant timescales as computed according to Eqs. (33)-(37) are also shown 
(right panel). Calculations were performed using the SkyNet nuclear network [39] (Courtesy of D. 
Vescovi). 
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where the sums run over all possible reactions that include (A;, Z;) in the initial or in 
the final state: the first sum contains decays, photodisintegrations and semi-leptonic 
processes, as electron, positron or neutrino captures; the second and the third ones 
include nuclear fusions with two and three reactants, respectively. In this context, a 
three-body reaction is a sequence of two-body reactions with an intermediate state 
with an extremely short lifetime. The factors ./' account for multiplicity effects in 
the case of identical particles: if Nm represents the number of m nuclei involved in 
a specific reaction with sign (i.e., Nm > 0 for creation and Nm < 0 for destruction), 
then Ni =N;, Nir = N;/(|N;|!|N;|!) and Nie =N;,/(|Nj|!|Ne|![N7|!) (for identical 
reactants, double counting must be avoided so that N; +N, =2 and Nj; +N; +N; =3 
for two- and three-body reactions, respectively. For example, for a-+a@+a— !C, 
N; =3 and N; = N; = 0). The calculation of the reaction rates and of the effective 
decay constants require the knowledge of the evolution of the fluid density and tem- 
perature, as well as information about the neutrino irradiation fluxes. In actual com- 
putations, the evolution of the matter density is usually prescribed or extracted from 
hydrodynamics simulations. The subsequent evolution of the temperature and Y, are 
then self-consistently determined from the detailed abundances, from the EOS of the 
plasma, and assuming that the expansion proceeds adiabatically, unless for nuclear 
energy generation and neutrino leakage or irradiation. 

We will now describe in more details the most relevant features of the nucleosyn- 
thesis happening in a fluid element expanding after a compact binary mergers, ac- 
cording to the outcome of nuclear network calculations. For the bulk of the neutron- 
rich ejecta, characterized by low specific entropy, we can identify four phases: 
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Fig. 5 For the same trajectory used for Figure 4, we also present detailed nuclear compositions 
on the nuclear chart. In this two panels we show the abundances at the end of the NSE phase (top 
panel), and the end of the (n, y)-(y,n)) equilibrium (bottom panel). In the latter case, the r-process 
path is clearly visible. These pictures were produced using the SkyNet nuclear network [39] and 
the dedicated visualization software (Courtesy of D. Vescovi). 


e the initial NSE phase; 

e the r-process nucleosynthesis phase; 
e the neutron freeze-out phase; 

e the decay phase. 


The different phases and their properties can be observed in Figures 4, 5 and 6. 


r-process nucleosynthesis from compact binary mergers 21 


2 22 Temperature = 3.00E+08 K 2 
5 bis aa » Density = 4.17E+01 g/cm? w 
— >k 08 
4 pig/cm] 1s al 
— &na ferg/s/ 9] Heating rate = 1.79E+17 erg /5/ 9 
o 16 100 


Entropy = 8.28E+01 KB / baryon „ ® 


é 
3 
c é 
S 2 Ye=0.300 ʻ 
Ed 2 
= 3 
z bg Abundance 
r T92 
É 
s 
3 a 
3 10 
= 1 
eo Unstable nuclide 107 
352 
a 0 Stable nuclide 
A A 107 
= Missing nuclide 
11 Closed neutron shell 10? 
= Closed proton shell 
Time = 1.16E+00 s 107 
=125 
10” 
10? 
10% g 10° 
f~™ 10° & 
3 ye 
| | w E 10 
as H J 3 
40 10 
36 02 10 
z T ë 
o 10? 
Made with SkyNet by Jonas Lippuner 2 72 94 111, 145 170 188 206 253 337 
Adapted by Diego Vescovi neutrons Ipeak  2™ peak rare-carth 3" peak fission material A 10? 


32 22 Temperature = 2.83E+05 K 
È g — TK 20 i 3 
= a = pig/m a Density = 6.54E-26 g/cm 
E — Sra ferg/sial |5 Heating rate = 7.43E+05 erg /5/ 9 102 
— = 100 
Se a Missed 3 Entropy =8.08€+02kB/baryon „ %® x 
` ug 6 220 
a -8 12 9 Ye=0.394 = ai 216 
a s 
2 -12 10 è 204 
~ -16| 200 
S 8 Abundance 
R -20 am 10? 
' T92 
É -24 188 
s 16 
S -28 19 E 
3 103 io? 10 10> 10° 10? 
œ% Time [s] 
aso 164 . a 
ase Unstable nuclide 10 
352 
ias 0 Stable nuclide 
-5 
i à 10 
= Missing nuclide 
5 
11 Closed neutron shell 6 
Shpucketors/ilippenerisrynet 10 
= Closed proton shell 
Time = 1.00E+09 s 107 
=32yr 
10” 
88 10? 
z s Fos 10“ 8 10° 
pee 10° 5 z 
u 0.4 | NESAS 10° E 10 
H A = 
502 gns 
! All | ae 
s ig o 7 10 
2 Made with SkyNet by Jonas Lippuner 2 72 94 111, 145 170 188 206 253 337 
ons +N ‘Adapted by Diego Vescov! neutrons 1 peak 2 peak rare-earth 3" peak fission material A 10? 


Fig. 6 Same as in Figure 5, but during the decay phase, just after neutron freeze-out (top panel) 
and at 10°s after merger (bottom panel). These pictures were produced using the SkyNet nuclear 
network [39] and the dedicated visualization software (Courtesy of D. Vescovi). 


The NSE phase 


The NSE phase (top panel of Figure 5) lasts as long as the temperature in the plasma 
stays above T ~ 4-5 GK. It is important to notice that a fluid elements can enter 
and exit the NSE phase several times before eventually decreasing its temperature 
below the NSE threshold due to hydrodynamics processes and intense nuclear heat- 
ing. In this case, only the conditions at the last NSE freeze-out influence the sub- 
sequent evolution. During NSE, since (p,y)-(y,p) and (n, y)-(y,n) reactions are 
all in equilibrium, the recursive application of the corresponding equilibrium rela- 
tion among the relativistic chemical potentials, namely U4 z) = Hp + H(a—1,z-1) and 
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Haz) = Un + U(a-1,z), Yields to U4 z) = ZUp +N Un. For particles described by the 
Maxwell-Boltzmann statistics, once p, T, and the abundances of free n’s and p’s 
(Yn,p), are provided, the abundance Y(,4 z) of any nucleus (A,Z) in NSE is given by 
[10]: 


3/2 A-1 3(A—1)/2 
Yon ayzZy 4-2) Guz) (T)A? / p 27h? eBIAZ)IkeT (24) 
(AZ) = pon 2A mp mpkgT i 


where B(4 z) is the nucleus binding energy. Yn,» are ultimately set by requir- 
ing baryon conservation and charge neutrality, ie. 1 = }(4 z) Y(a,z)A and Ye = 
Laz) Y(a zZ respectively. Very high densities and not too large temperatures fa- 
vor large nuclei (as it happens in the crust of cold NSs), while photodestruction in 
hot environments produces light nuclei and ultimately free protons and neutrons (as 
in the mantle above proto-neutron stars in core-collapse supernovae (CCSNe) or in 
the cores of merging NSs). For intermediate regimes, the nuclear binding energy 
favors the most tightly bound nuclei, i.e. iron group nuclei or œ particles among the 
light nuclei, and the temperature regulates the width of the distribution. 

For a given plasma configuration (i.e., for a given set of p, T and Y,) the NSE 
condition determines the nuclear abundances according to Eq. (24) without the need 
of solving Eqs. (23). However, the abundances change as a function of time due to 
the temporal evolution of the expanding and cooling plasma, p = p(t) and T =T (t). 
Moreover, weak interactions are out of equilibrium and Eq. (23) and the charge 
neutrality conditions reduce to an equation for Y,: 


dY, 
T = L (Aet i = Ae- i +h,.i — Aiei +Ag- — Ap+ i) Yis (25) 


l 


where the sum runs over all nuclei and the rates span all possible (if any) semi- 
leptonic reactions involving each nucleus i. 

As visible in the top panel of Figure 5, at NSE freeze-out the composition of the 
expanding plasma is characterized by a distribution of seed nuclei, peaking around 
the iron group, & particles and free neutrons. Depending mainly on the initial Y, < 
0.4, the seed nuclei are possibly very or extremely neutron-rich, with A ~ 60 — 100. 
A high fraction of free neutron is also expected. 

If the fluid is radiation dominated, its specific entropy can be approximated by 
the entropy of the photon gas, see Eq. (7). It is very insightful to substitute this 
expression inside the expression of the abundances in NSE conditions, Eq. (24), to 
express it also in terms of the entropy: 


3/2 3(A-1)/2 
Zy(A-z)_A kgT B(A,Z)/kgT 
Yaz) = Y Yn AA] (= eP(4Z)/kBT (26) 


Since Y4 z) x s-(4-1), we recover the result that in the case of matter with high 
specific entropy an a@-rich freeze-out is obtained. 
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The r-process nucleosynthesis phase 


Once T < 4 GK, the first reactions that run out of equilibrium are charged nuclear 
reactions involving the less abundant nuclei. Under the assumption that neutron cap- 
tures, photodestructions and B~ decays (possibly emitting j delayed neutrons, with 
j = 0 being the classical B~ decay) are the most relevant reactions, Eq.(23) be- 
comes: 


d¥(4,z) 
dt 


~ Mn (OV) (41,2) (n,y)(A,Z)(A—1,Z) FAA Zna z Yaz) + 


E [nn (ov) (A,Z)(n,y)(A41,Z) + NaAz)(yn)(A- 1 z] Yaz) T 


J 
+È Aasjz—y(az)te-tietin Y(a+j,z-1) + 
j= 


—V Ma2- Z+1)+e-+0+j n Yaz): (27) 
j=0 


In the first line, we have considered the creation of (A,Z) nuclei through (n, y) re- 
actions on (A — 1,Z) nuclei and (7,7) reactions on (A + 1,Z) nuclei. In the second 
line, we have considered the destruction of (A, Z) nuclei through (n, y) and (7,7) re- 
actions. In the third and forth lines, we have taken into account B~ decays that can 
include a (A,Z) nucleus in the final or initial state, respectively, possibly through 
the additional delayed emission of 0 < j < J neutrons. 

As we have seen in the previous section, if the temperature and density are large 
enough, the timescale of B~ decays are much longer than the timescales of (n, 7) 
and (y,n) reactions. The latter absorption processes connect nuclei among the same 
isotopic chain (Z is fixed) and their equilibrium ensures an almost steady free neu- 
tron fraction, while the former decays connect nuclei of contiguous isotopic chains 
(Z— 1 and Z, and Z to Z + 1). Then, the evolution predicted by Eq. (27) can be split 
into two separate problems: the much faster (n, y)-(n, y) equilibrium inside each iso- 
topic chain; the slower flow through different isotopic chains driven by p~ decays. 
In the following we will closely analyse each of the two problems separately. 

We start by consider an isotopic chain characterized by a specific Z. The equilib- 
rium condition inside the chain (i.e. H(4+1,z) = U(a,z) + Hn) allows to write an equa- 
tion for the abundances of two adjacent nuclei in the chain, starting from Eq.(18) 
and simply noticing that Y(441,z)/Y(a.z) = n(4+1,z)/n(a z): Evaluated the resulting 
expression for typical magnitudes, we obtain: 


Y(a+iz) ~ 571x103 Nn iat \ OP exp (S,/kpT ) 
Yaz) 8 x 1074cm~3 J \ 0.1MeV exp (30) 


We notice that the abundance ratio depends only on nn, T and on S,,. The latter intro- 
duces a dependence on the nuclear masses. Also in this case, the numerical value of 
this ratio is dominated by the exponential factor whose argument compares the neu- 
tron separation energy with the plasma temperature. Close to the valley of stability 
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Sn ~ 15MeV > kgT and Y(441,z) > Y(4,z), i.e. abundances increase steeply moving 
toward neutron richer nuclei. Assuming kgT = 0.1 MeV and S,, = 3MeV, the above 
estimate still implies Y/4,1z) >> Y(4,z). For Sn = 2MeV, the numerical prefactor de- 
creases to 0.26, meaning that Y(441,z) < Y(4,z). Since Sn — 0 at the neutron drip 
line, there exists always a turning point in the isotopic chain, where Y(4+1 z) stops to 
increase with respect to Y(4 z) before decreasing when approaching the neutron drip 
line. We estimate this point by requiring ¥(441,z)/Y(4,z) = 1. Assuming a certain nn 
and temperature T, this condition translates in a reference neutron separation energy 


50: 
2 /mpgkgBT\?” kpT 
0 Bkp B 
= keT in | — = 2.14 Mev ( “B= 
Sa = k a(ż ( nh? ) ) e (w) i 


Nn kgT 


The value of s? depends on n, and T, but not on Z. So, at any given time, all pop- 
ulated isotopic chains have their abundance peaks at nuclei characterized by the 
same $? and the conditions S,(A,Z) = S? defines the r-process path. A better ap- 
proximation can be done by considering that, due to nuclear paring effect, even 
neutron numbers are favored. Indeed, while S, decreases for increasing D, but with 
a even-odd modulation, S2,, defined as $3,(A +2,Z) = (m4 z + 2m, — m4 42,2))C” 
and called the two neutron separation energy, decreases in a smoother way, with a 
sudden decrease at the magic neutron numbers. By considering the effective equi- 
librium (A + 2,Z) +y 4+ (A,Z) +2n and the corresponding relation on the chemical 
potentials, U442,7) = H(a,z) +2Hn, the r-process path can be defined as the set of nu- 


clei belonging to different isotopic chains for which S$2,(A,Z) is closer to 25° (ny, T) 
along each specific chain. 

Since the relative distribution of abundances inside an isotopic chain is set by 
the fast (n, y)-(y,n) equilibrium and determined by Eq. (28), it is useful to consider 
the total abundance along the chain, Yz = )’4 Y(4,z). Starting from Eq. (27) and con- 
sidering only B~ decays without delayed neutron emission (J = 0), we obtain an 
evolution equation for Yz: 


dýz > x 
es Ap (z—1)¥z-1 — Àg zYz (30) 
where Ap.z’s are the effective, abundance-weighted B decay rates of the whole 
chains: 


Agz = [Easieaias] / (Eran) . (31) 
A A 


Under the assumption that the amount of free neutrons and the temperature are 
high enough, the duration of the r-process becomes comparable or even larger than 
the B decay lifetime (Eq. 20) over the most relevant part of the nuclear chart, and 
the chain abundances Yz tend toward an equilibrium configuration, i.e. dYz/dt ~ 0, 
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which translates into: : : 
AzYz = Àz- 1Yz—ı X const. (32) 


This condition is knows as steady B -flow approximation [33, 16] and it implies that 
the abundance of a chain is proportional to its effective B decay lifetime. Since 
the nuclei with magic neutron numbers N = 50,82,126 (or just above them) and 
closer to the valley of stability have the longest B decay lifetime, it is expected that 
maxima in the abundances will occur at the top end of the kinks in the r-process 
path corresponding to the neutron shell closures. This is what quantitatively defines 
the location and the relative importance of the waiting points inside the r-process 
path. 

The duration of the r-process nucleosynthesis phase crucially depends on the 
availability of free neutrons. Under the assumption that neutron captures and pho- 
todestructions are the most relevant reactions, the timescale over which the neutron 
abundance changes, T,, can be evaluated from Eq.27 as: 


1 1 1 
; (33) 
Yn/Ysced (2 ) Ty,n) 


The relevant timescales T, y) and Tyn) can be estimated through Eqs (15) and (19) 
as we did in the previous section, but they can also be computed more rigorously 
from the reaction rates that enter Eq.(23): 


1 1 
Tn Y, 


dY, 
dt 


1 
“(5 Yaz Mn (OV) (A.Z) (ny) wasn) Han (34) 
( 


Tn,7) AZ) 


1 
=( Z rza \(y,n)(4— a)/) È Yaz: (35) 
( 


Uyn) AZ) 


As long as the neutron-to-seed ratio is large, T, is relatively long and there is enough 
time for many isotopic chain (with increasing Z) to be reached and for the B-flow 
equilibrium to establish. Detailed calculation show that, for typical conditions in 
compact binary merger ejecta, this r-process phase can last up to ~1 second af- 
ter the merger (see for example the bottom panel of Figure 5). Eventually, the r- 
process path can reach the neutron magic number N = 184. Above that, fission 
becomes the dominant nuclear process. Typical super-heavy, neutron-rich nuclei 
have (A,Z) ~ (250,100), while daughter nuclei can be approximated by distribu- 
tion around the double-magic number nucleus (Aj 2,Z1,2) ~ (132,50). Thanks to 
fission and depending on the lighter fission fragment distribution, the abundances of 
heavy elements double and specific region of the nuclear chart start to populate. If 
the initial neutron richness is high enough (in particular, when the initial Y, < 0.1) 
fission becomes relevant while (n, y) and (7,7) reactions are still in equilibrium. The 
lighter nuclei resulting from the fission still capture many neutrons, reaching again 
N = 184 and many fission cycles occur [32, 43, 38]. 
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The possibility of reaching (n,y)-(y,n) equilibrium and its duration are criti- 
cally related with the temperature evolution. For a radiation dominated plasma, 
if the expansion proceeded adiabatically the temperature would simply evolve as 
T(t) © Thom(thom/t), with Thom of the order of a few GK. Then, within ~ 0.1 s, 
the temperature would drop below 1GK. While this relation describes in good ap- 
proximation the evolution of matter temperature before and after the r-process, it 
is rather inaccurate during it. Indeed, during the r-process the energy released by 
nuclear reactions and the matter density are often large enough such that matter 
is significantly re-heated and temperature becomes again larger than 1 GK for a 
much longer time (up to ~ 1 s), before dropping again as predicted by the adiabatic 
expansion law, T œ t~l. To estimate if the specific nuclear heating épuc) (nuclear 
energy per unit mass, per unit time) can affect significantly the matter temperature, 
we compare the density of energy of the radiation field with the energy released 
by the nuclear decays during a time comparable to the expansion timescale 7, i.e. 
nkt Tax /15(hc)? ~ énuciP T, and we estimate énuct ~ SnYn/(mpAtr—proc) [43]. By 
solving for Tmax, We can estimate the maximum temperature reached by matter due 
to the nuclear heating. Assuming that around t ~ 0.1s the density has decreased 
down to 10°g cm~? and the (n, y)-(y,7) equilibrium lasts for At,—proc ~ 1 8: 


Y, 1/4 ] 1/4 1/4 (At, -1/4 
Tmax = 0.75 GK P a (—_) ey 
0.5 x 10° g cm~? 5 MeV 10 ms ls 
(36) 


When Tmax Z 0.7 GK the temperature is large enough to guarantee (n, y)-(Y,n) equi- 
librium during most of the neutron captures. This is defined as the hot r-process. If, 
on the other hand, Tmax < 0.3 GK, the density is large enough to sustain the neutron 
captures, but not (n, y)-(y,n) equilibrium. In this latter case, the neutron density 
decreases much faster due to neutron consumption and this is referred as cold r- 
process. It is important to stress again that the relevant r-process regime (hot or 
cold) does not depend on the absolute peak temperature, but on the maximum tem- 
perature during the re-heating phase occurring when seed nuclei capture neutrons. 
Other factors, as for example the expansion timescale, determined if the r-process 
happens in hot or cold conditions. 


The neutron freeze-out and the decay phases 


As the r-process proceeds, and matter expands and cools below ~1GK, Tyn) > 
Tiny) and photodestruction becomes inefficient in keeping a high neutron density 
outside the heavy nuclei, while the still effective neutron captures and B~ decays 
produce more and more heavy nuclei. Once Y,,/Yseea ~ 1, according to Eqs. (17) 
and (33), it happens that T, ~ Y,,/ Yseed T(n,y) < 1s: the free neutron lifetime suddenly 
decreases because the seed nuclei start to compete for the few available neutrons. 
This phase is called neutron freeze-out and it is characterized by a sudden drop of 
the neutron density and of the neutron-to-seed ratio, as visible both and in Figure 4 
and in the the top panel of Figure 6. After this drop, the B~ decays (often followed 
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by the emission of free neutrons) start to compete with the neutron capture, since 
Tín, y) becomes comparable to Tg. The latter can be estimated as Eq.(20), but also as: 


1 
—= ( y y Yazihasacizane sei] i( y Yaz) : 
(A,Z) (A 


,A>4 j=l J Z), A>4 
(37) 


The presence of freshly emitted free neutrons provides a new source of neutrons 
available to be captured also at later time. The result is that, while the nuclei that 
were located along the r-process path during the (n, y)-(y,n) equilibrium decay col- 
lectively toward the valley of stability, the competition between the neutron captures 
and the B decays smooths the r-process abundances. In particular, it removes the 
strong oscillations in the mass number that characterize the abundances just before 
neutron freeze-out, due pairing and collective effects in the nuclear properties. 

After a timescale ranging from a few up to a few tens of seconds (depending 
on the environment and increasing for decreasing initial Y.) the neutron density has 
decreased such that also the neutron captures become negligible. Most of the matter 
is still in the form of unstable nuclei, with neutron excess D ~ 5-10. At this stage 
B~ decays and, depending on how extended the r-process path is, œ decays are 
the most relevant nuclear processes that bring the abundances toward the valley of 
stability. While most of the nuclei have reached stable configurations within a few 
tens of days, a few heavy isotopes have longer half-lives, extending above 10” yrs 
for 14 isotopes and above 10° yrs for 7 isotopes. Since all these reactions liberate 
nuclear energy, they can heat up matter or become a distinct source of radiation. The 
detection of the various, distinct signals is thus a clear indication of the r-process 
nucleosynthesis occurring inside the ejecta. 


The r-process peaks and the s-process nucleosynthesis 


During the (n, y)-(y,n) equilibrium the r-process path is characterized by the pres- 
ence of waiting points, where matter accumulates due to the longer lifetimes oc- 
curring in correspondence of the neutron magic numbers N = 50, 82,126. Waiting 
points are identified by neutron close shells, but they can span a large range of 
Z’s and they can also be influenced by the presence of the proton magic numbers 
Z = 28,50 (see the bottom panel of Figure 5). Since at neutron freeze out nuclei 
are still very neutron-rich (typically Z/A ~ 0.5), these waiting point nuclei are char- 
acterized by A ~ 80, 130, 195. When decaying back to stability, nuclear abundances 
proceed along A ~ const paths, i.e. they keep their mass number approximately con- 
stant while increasing their Z. However, the residual delayed neutron emissions and 
neutron captures can further change A by a few units during the process. Due to the 
decreasing neutron density, the effect due to B decays eventually takes over. The 
result is that the final abundance pattern is characterized by three peaks (known as 
the r-process peaks): the first peak is located around A ~ 80 and covers the range 
A ~ 72 — 94, the second around A ~ 130 with range A ~ 111 — 145, and the third 
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around A ~ 195 with range A ~ 188 — 206. Due to the B decays, the nuclei neutron 
content has globally decreased, the Y, has increased toward < 0.4. As a function of 
Z, the peaks are located around Z ~ 35 (e.g. Selenium, Bromine, Krypton) for the 
first peak, Z ~ 53 (e.g. Tellurium, Iodine, Xenon) for the second peak, and Z ~ 77 
(e.g. Osmium, Iridium, Platinum) for the third peak. Thus, due to the combination 
of multiple neutron captures followed by many f decays, the final peaks are shifted 
toward the left by several units with respect to the nuclei characterized by magic 
neutron numbers along the valley of stability. 

However, when considering the element above Zn, the abundances that we ob- 
serve inside the solar system reveal, close to the r-process peaks, other peaks shifted 
toward the right by a few units in the mass number and happening precisely in cor- 
respondence of the magic neutron numbers. Then, differently from what happens 
in the r-process, the production of these peaks must make place close to the valley 
of stability, in such a way that a neutron capture on a stable or long-lived nucleus 
is followed by a B decay before another neutron can be captured, meaning that 
t(n,y) > tg. This is possible if the neutron densities are in the range 10° en’, 
so several orders of magnitudes lower than the ones required by the r-process. This 
kind of nucleosynthesis is called slow neutron capture process (s-process) and it is 
thought to happen inside low and intermediate mass stars (starting their life with a 
mass between 0.6 and 10 Mọ) during their asymptotic giant branch (AGB) phase. 
Starting from iron group nuclei produced in a previous SN explosion and already 
present inside the star, the s-process proceeds through the s-process path up to lead 
on a timescale of several thousands of years. The abundances of the s-process nu- 
cleosynthesis are well understood since they rely on the well-known properties of 
nuclei at or very close to the valley of stability. Thus, the solar r-process abundances 
are defined as the residual of the solar abundances, once the s-process contribution 
has been removed, see e.g. [51]. Most of the heavy elements receive a contribution 
from both the s- and the r-process. It is important here to stress that information 
about the isotopic composition is in general harder to obtain, while astrophysical 
observations often provide information on the elemental composition. So, an ele- 
ment usually receives contributions from different isotopes and these isotopes can 
be synthesised in different environments. Nevertheless the structure of the valley of 
stability is such that some stable isotopes can only be produced in the s-process or 
in the r-process. In the first case, this happens when a stable isotopes is shielded 
by another stable isotopes along a B decay line. In the second case, this happens 
when a stable nucleus can be reached by a neutron capture and a subsequent B de- 
cay sequence only from a nucleus that has a very short lifetime or that is outside 
the s-process path. When this happens for all or for the most relevant isotopes of 
an element, this element is called s-only or r-only element, respectively. A r-only 
isotope has the advantage of allowing the study the contribution of the r-process 
independently from possible contamination from the s-process. Europium and el- 
ements above lead (including Uranium, Thorium and Plutonium) are among the 
r-only elements. 
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Nucleosynthesis in high entropy and fast expanding ejecta 


Not all the ejecta are characterized by low or moderate entropy, or by expansion 
timescales ranging from a few ms up to several tens of ms. In particular, detailed 
merger models show the presence of high entropy and/or extremely fast expanding 
tails in the ejecta distribution. Moreover, it is not always guaranteed that the peak 
temperature will be above 4GK or that the temperature will be high enough during 
the r-process phase to ensure (n, y)-(y,n) equilibrium. The nucleosynthesis in these 
conditions can be significantly different than the one described above. 

For example, if the expansion timescale t becomes of the order of (or even 
smaller than) 1 ms, the re-heating during the r-process phase could become inef- 
ficient and (n, y)-(Y,7) equilibrium cannot establish, see Eq.(36). In this case, a cold 
r-process happens. Moreover, despite the possible large abundance of free neutrons 
and the availability of seed nuclei, the decrease in the neutron density is so fast that 
neutron captures become inefficient too early and the r-process does not proceed 
much. A large fraction of the free neutrons are not captured by nuclei and decay 
into protons. The resulting abundances are very different from the one observed in 
a full r-process. According to Eq.(5), for the dynamical ejecta emerging from BNS 
mergers, this is the case for matter expanding at v ~ 0.6 — 0.8c. Such ejecta are ob- 
served in some models in the high speed tail, especially if the EOS of nuclear matter 
is rather soft and the shocks produced by the NS collision are very violent. 

Finally, if the specific entropy is large enough (s > 80 kgbaryon~ '), the NSE ends 
with a a@-rich freeze-out (see e.g. [29] and reference therein). If, additionally, Ye < 
0.4, there is a significant abundance of free neutrons that are not bound inside the 
very abundant & particles and inside the fewer seed nuclei. Thus, in these conditions 
the neutron-to-seed ratio is still very high and the r-process nucleosynthesis can still 
occur to produce heavy r-process nuclei. The fundamental different is that the final 
abundances are overall dominated by the œ particles emerging from the NSE freeze- 
out, see e.g. [52, 21]. 


Nuclear physics input and detailed network calculations 


Equilibrium argument are very useful to get a qualitative understanding of the main 
features of the r-process nucleosynthesis. However, only the numerical solution of 
Eqs.(23) provides detailed nucleosynthesis yields. In nuclear network calculations 
abundances are typically initialized in NSE conditions at high temperatures (T ~ 6 — 
8GK) and then evolved consistently following the chemical evolution of the ejecta 
during their expansion. For r-process nucleosynthesis calculations more than 7000 
nuclei are necessary, ranging from free neutrons and protons up to very neutron-rich 
transuranic elements (e.g. Curium), and covering the whole neutron-rich side of the 
valley of stability. The reaction rates needed by the network require a vast amount 
of nuclear physics knowledge. In the following, we briefly discuss the most relevant 
inputs, namely the nuclear masses, the B-decay half lives, the neutron captures rates 
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and the fission physics. For more detailed information we address the reader to 
Ref [11]. 

Nuclear masses. The most basic nuclear property are the masses of all nuclei, 
since they determines the threshold energies for all relevant reactions and, in par- 
ticular, the neutron separation energies. While the masses of stable nuclei and of 
nuclei close to the valley of stability are experimentally well known, only theoret- 
ical values predicted by nuclear mass models are available for exotic neutron-rich 
nuclei. Nuclear mass models are tuned on experimentally known masses and then 
extrapolated, so that the uncertainties grow moving toward the neutron drip line and 
just above the shell closure points, where correlations and deformation may be very 
relevant, but difficult to be taken into account. 

Beta decay rates. Also the values of p decay half-lives are experimentally un- 
known far from the valley of stability and theoretical calculations are necessary. 
Their values are crucial, since they determine the matter flow between different iso- 
topic chains. In particular, the most relevant decay rates are those of nuclei at the 
neutron magic numbers, since they are the waiting points of the r-process path. The 
calculation of these decay rates dependes on the reaction Q-value (and thus on the 
mass model), on the transition strength and, in particular, on its energy distribution. 
If the transition leads to a final state whose energy is above the neutron separation 
energy, the emission of one or more neutrons in the final state (B-delayed neutron 
emission) is possible. This emission is relevant after the neutron freeze-out phase. 
Also in this case, shell effects and nucleon correlations are possibly very relevant 
but difficult to model. 

Neutron capture rates. While neutron captures are essential to establish (n, y)- 
(y,n) equilibrium, the corresponding QSE condition makes the detailed knowledge 
of their rates not very relevant in modeling the r-process. It is after neutron freeze- 
out, when the temperature drops below | GK and the photodisintegration becomes 
less and less relevant, that a detailed knowledge of the relevant neutron capture rates 
becomes more important. At this stage, nuclei are still very exotic and theoretical 
calculations are once again necessary. They are usually computed within the sta- 
tistical model of nuclei, which is done using the Hauser-Feshbach approach and 
modeling the nuclear energy density, the y-strength function for the decay of the 
compound nucleus, and the potentials of several light particles. 

Nuclear fission. The fission properties of super-heavy nuclei (A < 280) are 
very relevant in calculations involving very neutron-rich matter (at least for ini- 
tial Y, < 0.10), but they are very uncertain. The theoretical description of fission is 
highly non-trivial, and it mainly depends on the fission barriers, defined as the en- 
ergy required for a nucleus to undergo fission. This energy is necessary to deform 
the nucleus such that the transition to the fission fragments states becomes ener- 
getically favored. If the transition to the final state does not require any additional 
energy, the fission is spontaneous. If additional energy is provided by the interaction 
or the absorption of another particle, then the fission is induced. All possible fission 
channels (e.g., spontaneous, neutron-induced, y-induced) and all the other compet- 
ing reactions (e.g. @ and p decays) need to be included in the model to provide 
reliable predictions. In the case of a very neutron-rich environment, as the one typ- 
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Fig. 7 Abundances as a function of the mass number (A, left panel) and of the atomic number (Z, 
right panel) at 10° seconds after merger for trajectories characterized by s ~ 11kgbaryon™! and 

~ 11 ms, but for different initial Y,’s, computed using the SkyNet nuclear network [39]. Black 
dots represent the Solar r-process residual, as reported by [51] (Courtesy of D. Vescovi). 


ically expected in the ejecta of compact binary mergers, neutron-induced fission is 
usually the most relevant channel. The very uncertain fission fragment distributions 
are also very important as they influence the abundances around A < 140. 


While the basic features of the r-process nucleosynthesis in the ejecta of compact 
binary mergers are robust and, especially for low entropy ejecta, mainly influenced 
by the initial Y,, the nuclear input physics can significantly change the fine struc- 
ture of the abundance pattern. For example, the B decay half-lives of nuclei with 
Z = 80 regulate the mass flow toward the magic neutron number N = 126 (possibly 
enhancing the amount of nuclei that go though fission cycles) and possibly affect 
the position of the third r-process peak [13, 43]. Moreover, since fission fragments 
are located around A < 140, fission physics can influence the width of the second 
r-process peak, as well as the abundances of the rear-earth elements, between Lan- 
thanium and Lutetium [19]. 


Detailed network calculations and nucleosynthesis yields from 
compact binary mergers 


Once all the relevant nuclear input physics has been considered, nuclear network 
calculations can predict with accuracy the distribution of the yields synthesised by a 
specific fluid elements. Nuclear network calculations for fluid elements expanding 
adiabatically show that the nucleosynthesis outcome depends mainly on the three 
parameters that we have used to described the properties of the ejecta at the onset of 
the homologous expansion phase, namely the electron fraction, the specific entropy 
and the expansion timescale (see, e.g., [38] for a systematic study). As representative 
cases for the abundances obtained at ~ 30y after merger, in Figures (7) we present 
the abundance patterns obtained by detailed nuclear network calculations. All trajec- 
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Fig. 8 Mass integrated abundances in the dynamical ejecta obtained by simulations of BNS merger 
models. In the left panel, the different curves represent different neutrino treatment,with the blue 
curve not including neutrino irradiation, while the other two including it. In the right panel, neutrino 
irradiation is always accounted for, but the different curves show the difference between an equal 
(blue) and an unequal (yellow) mass mergers. Figures taken from [54]. 


tories have the same specific entropy and expansion timescales (s ~ 1 1kgbaryon7! 


and T + 11 ms), but differ because of their initial electron fraction. In all cases, 
we compare the calculation outcome with the solar r-process residual. For Y, < 0.2 
all elements between the second and the third r-process peaks are synthesized with 
a pattern that well reproduces the solar one. Here the observed differences in the 
shape of the rare-earth peak could be mostly ascribed to the nuclear physics input. 
This nucleosynthesis outcome is often referred as the ’strong r-process”. For very 
low electron fractions (Y, < 10) the nucleosynthesis proceeds up to Ur and Th, and 
the abundance pattern is very robust, due to the extremely long r-process path and 
to the effect of several fission cycles, while for 0.1 < Y, < 0.2 the production of 
actinides decreases. Around Y; crit 0.23 — 0.24, the neutron-to-seed ratio at NSE 
freeze out is only of a few, the production of nuclei above the second r-process peak 
is inhibited, and the abundances move toward the first r-process peak for increasing 
Yo. This nucleosynthesis outcome is sometimes called ’weak r-process”’. 


Nucleosynthesis yields from compact binary mergers 


Detailed abundances in the ejecta from compact binary mergers reflect the ranges 
and the relative relevance of the distributions of entropy, expansion timescale and 
electron fraction that emerge from the merger dynamics. Different abundance pat- 
terns characterize the different ejection channels and keep an imprint of the merger 
nature, dynamics and aftermath. 

Dynamical ejecta. The low-entropy, equatorial dynamical ejecta have Y, < 0.2 
and produce robustly all elements between the second and the third r-process peak. 
In addition, if Y, < 0.1, significant abundances of translead nuclei and actanides 
are produced in relative abundances that can be comparable or even larger than the 
solar ones. This is the case especially for BHNS mergers and BNS mergers char- 
acterized by very different NS masses. In the case of BNS mergers of comparable 


NS masses, at polar latitudes the entropy increases but stays on average below 20 
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Fig. 9 Mass integrated yields in the dynamical ejecta obtained by simulations of BHNS merger. 
In the left panel, mass fractions are presented for different combinations of NS and BH masses. In 
the right panel, the different curves represent abundances obtained for different intensities of the 
neutrino luminosities. Figures taken from [29] (left panel) and [55] (right panel). 


kgbaryon™!, while Y, > 0.25 and the production of heavy r-process elements is sup- 
pressed. In this case, the r-process nucleosynthesis does not proceed up to very high 
mass numbers, but produces elements from the first up to the beginning of the sec- 
ond peak. The relative importance of the two contributions (polar and equatorial) to 
the dynamical ejecta depends on the binary mass ratio and on the properties of the 
nuclear EOS. However, the equatorial component is expected to be always present 
and overall dominant. Representative results for mass integrated abundances in the 
dynamical ejecta of BNS and BHNS simulations are presented in Figures 8 and 9 
and they show that the dynamical ejecta can produce a significant fraction of all 
r-process nuclei from the first to the third peaks, see e.g. [70, 54]. 

Disk wind ejecta from BH-torus systems. For a BHNS merger or for a BNS 
merger whose the central remnant collapses to a BH within a few dynamical 
timescales, the ejecta from disk winds are dominated by the viscous component, 
have low entropy and cover a broad range of electron fractions across the critical 
value Y; crit ~ 0.23 — 0.25. In these conditions, the production of all heavy elements 
between the first and the third r-process peak is expected (see e.g. [29, 76]). Due 
to the smaller effect of neutrino irradiation, the synthesis of elements between the 
second and the third r-process peaks is significant while the production of the first 
peak elements is below the solar ratio. The angular distribution of the abundances 
is in this case rather insensitive to the latitude: high-Y, matter at polar latitude pos- 
sibly synthesize only elements below the second r-process peak, but this v-driven 
component is very subdominant for BH-torus systems, due to the lower neutrino 
luminosities. Nucleosynthesis results for the viscous ejecta from BH-torii systems 
are presented in the left panel of Figure 10. 

Disk wind ejecta from systems hosting a massive NS. If the massive NS re- 
sulting from the merger of two NSs is not a short-lived one, a larger variability is 
expected, reflecting the presence of several components in the disk winds and the 
impact of neutrino irradiation [37]. If the remnant survives for a timescale compa- 
rable to the disk lifetime (facc, see Eq.(3)), neutrino irradiation progressively shift 
the Y, distribution of the viscous ejecta toward larger values = Y% crit, compared to 
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Fig. 10 Mass integrated nucleosynthesis results (as a function of the mass number A) for the disk- 
wind ejecta obtained in the case of a BH-disk systems (left panel) and of a massive NS remnant 
(right panel, with different colors corresponding to different NS lifetimes). Figures taken from [76] 
(left panel) and [37] (right panel). 


the BH-torus system. While the production of both light and heavy r-process el- 
ements is forseen, the relative importance of the former with respect to the latter 
increases as a function of the remnant lifetime. The presence of a very long-lived 
massive NS (collapsing on a timescale > facc) could possibly prevent the formation 
of the second peak and of all the elements above it in the bulk of the ejecta, produc- 
ing a significant amount of light r-process elements. The presence of a massive NS 
causes also disk wind ejecta emerging from distinct portions of the solid angle to 
be characterized by different nucleosynthesis patterns. As in the case of the dynam- 
ical ejecta, matter expelled at high (polar) latitudes is more efficiently irradiated by 
neutrinos and for it the production of heavy r-process elements at and beyond the 
second r-process peak becomes soon inefficient, while elements characterized by 
75 SA < 120 and 33 Đ Z < 55 are eventually produced in solar proportions even 
for ejecta expelled within the first few tens of milliseconds [47, 40]. Nucleosynthesis 
yields obtained for the neutrino- and the viscosity-driven ejecta from merger rem- 
nant hosting a massive NS in the center are presented in the right panel of Figure 10. 


High entropy, high velocity tails. In the case of shock-heated ejecta and ejecta 
expelled from strong magnetic fields in torii around BHs, a high entropy tail is 
observed in the ejecta distribution. In the first case, the ejecta also expand with a high 
velocity, such that T < 2ms. In this case, assuming Y; < 0.4 (but often Y, < 0.3) the r- 
process nucleosynthesis proceeds in a@-rich freeze-out conditions and it synthesises 
a significant fraction of nuclei up to lantanides (Y, < 0.4) or even actinides (Yo S 


0.2). Differently from the low entropy, neutron rich ejecta, these matter present also 
a significant amount of H and He in the final composition. 
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Observables of compact binary merger nucleosynthesis 


Binary compact mergers have long been thought as promising astrophysical sites for 
r-process nucleosynthesis. While the expected conditions of the ejecta and the out- 
come of detailed numerical models indicate that the production of heavy r-process 
elements above the iron group is robust, observational evidences are crucial to vali- 
date our models and to discriminate between the many theoretical uncertainties that 
still affect our theoretical understanding. 

In the following we will discuss two major observables and their relation with 
the outcome of r-process nucleosynthesis in compact binary mergers: the electro- 
magnetic transient called kilonova and the evolution of the chemical abundances of 
r-process elements in the stars of our Galaxy and in its satellites. 


Electromagnetic signatures of r-process nucleosynthesis in 
compact binary mergers 


What is a kilonova? 


Starting from a few seconds and continuing for several hundreds of days after 
merger, a large amount of nuclear energy is released by the combination of p de- 
cays, œ decays and fission processes that follow the r-process nucleosynthesis in 
the ejecta of a compact binary merger. This energy heats up the expanding matter 
and produces a nuclear powered transient called kilonova [36]. Depending on the 
ejecta properties, a kilonova is expected to have its peak luminosity between a few 
hours and several days after merger in the UV/optical/near-IR frequences, with a 
fast declining luminosity. For a detailed treatment of this transient we refer to the 
dedicated Chapter and to recent reviews, e.g. [14, 45, 11]. 

The nuclear energy associated with each decay is distributed among the daughter 
nuclei and the other particles in the final state (electrons, photons and neutrinos), in 
a way that primarily depends on the nature of the decay. Typical Q-values are of the 
order of 1-100 MeV, much larger than the decreasing matter temperature in the ex- 
panding ejecta. Indeed, if the temperature at the end of the r-process nucleosynthesis 
(tr—proc ~ 1 8) is kg T;-—proc ~ 0.1MeV, then: 


t \! kpT,— t t \-l 
kgT (t) ~ keT;—proc (=) ~ 1.16 eV E <a) ( meo) (; -) 

(38) 
While the energy emitted in neutrinos is always lost, several processes can ther- 
malize at least a fraction of the released nuclear energy, making it available for the 
kilonova. The intensity of the nuclear heat and the efficiency of the thermalization 
processes are maximal immediately after the r-process nucleosynthesis. However 
matter is initially very opaque to photons and the timescale for photon diffusion is 


36 A. Perego, F.-K. Thielemann and G. Cescutti 


much larger than the dynamical timescale over which the ejecta expand. The pho- 
ton optical depth depends on the matter density profile and on the photon opacity. 
It is thus necessary to wait for the density to drop such that thermal photons can 
efficiently diffuse and be emitted at the photosphere. 


r-process nucleosynthesis and kilonovae 


If this is the mechanism behind kilonovae, r-process nucleosynthesis can influence 
these transients mainly through three aspects: the nuclear heating, the thermalization 
efficiency and the photon opacity. In the following we will analyze the main features 
of each of them. 

Nuclear heating. On the kilonova timescale nuclear abundances in the ejecta 
change in time according to a set of decays and fission reactions. Each reaction 
is characterized by an exponential behavior, N;(t) = Nioexp(—Ait) where N;(t) is 
the number of parent nuclei at time ft and A; the reaction rate. The specific nuclear 
heating rate as a function of time can be computed as: 


_ L; OiANi(t) 


one Z êe- (t) +éy(t) +éa(t) + éy(t) + fission (t) , (39) 
ej 


ênuci (t ) 
where Q; is the Q-value of each reaction i and mej is the ejecta mass. In the last step 
we have explicitly indicated the decay particles whose kinetic energy provide the 
available nuclear energy, neglecting the kinetic energy of the daughter nuclei in the 
case of œ and f decays. Detailed network calculations (e.g., [44, 32, 38]) show that 
the dominant contribution to nuci can be approximately described by a power-law 
dependence in time, with possible corrections in the form of exponential terms: 


N 
ênuci (t) X Enucl,0 (t/t) * +) Biexp( -y t). (40) 
i=1 


with fo being a reference timescale. The precise values of énucio, Œ, P; and y; depend 
on the initial ejecta properties, i.e. on (Y.,s,7) of the specific trajectory, and on 
the nuclear input physics. For low entropy ejecta, Y is the dominant parameter. 
For Y, < 0.25, œ ~ 1.3 with a possible variation interval 1.1 — 1.4, while different 
mass models give épucio ~ 10° — 10"erg s7! g`! (assuming tọ = 1 s) [57]. This 
is mainly due to the fact that, especially for very low Y,, translead nuclei tend to 
decay trough œ decays and fission, which are very sensitive to the specific mass 
model [6, 57, 75]. The dominant power-law term in Eq. (40) can be understood 
by considering that in most of the cases (and especially for Y, < 0.2, when strong 
r-process nucleosynthesis occurs) a large statistical ensemble of nuclei is produce 
[44]. As explained in [44], assuming B~ decay to be the dominant decay channel, 
along an isotopic chain (i.e. for fixed Z), the Q-value is roughly proportional to the 
neutron excess D. Since Ag œ Q°, then Ag œ D°. The nuclei at the end of the r- 


process distribute along the r-process path, characterized by S,, ~ S°, see Eq.(29). 


n? 
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According to the S, distribution on the nuclear chart, at high N a larger D is required 
to fulfill the r-process path condition and, for a given time f, the nuclear distribution 
is such that the number of nuclei per interval of D is approximately constant within 
the relevant interval of neutron excess: 


a Seon =CA—4/5 | (41) 


dA dDda 
with C being a constant. Due to the large number of reactions and nuclei involved, 
Eq.(39) can be converted in an integral over À: 


+0 
enai = [ AMA exp(—Arja T4 dA 1-7/9, (42) 


While this argument explains well the presence of the power-law and provides a 
good estimate of its slope, many details can affect the precise calculation, accounting 
for the lower exponent obtained in detailed calculations. 

The power law term in Eq.(40) is very robust for ejecta characterized by Y, < 0.25 
and undergoing strong r-process nucleosynthesis. For Y, = 0.25, the decay of a few 
nuclei dominate the heating rate and deviations from the power-low behavior are 
accounted by the exponential terms in Eq.(40). These nuclei (including for example 
neutron-rich isotopes of Kr, Rb, Br, Sr, Ce, As, and Ge, see e.g. [40]) are usually 
below the second r-process peak and are characterized by large Q-values and half- 
lives of a few hours. Thus the heating rate is enhanced by a factor of a few within 
the first day and decreased by a factor ~ 2 at later time with respect to the Y, < 
0.25 case. However, since the exponential corrections depend significantly on the 
ejecta distribution in the (Y.,s,7) space and the latter is usually rather broad in Y,, 
a power law behavior is partially recovered and nuc] varies only by a factor of a 
few for relevant ejecta conditions during the first week after merger (e.g. [38]). At 
late time (between 10 and several hundreds days), independently from the initial 
Y,, the nuclear distribution is very close to the valley of stability and only very few 
nuclei have the right lifetime to decay within this time window. Thus, the heating 
rate becomes possibly very sensitive to the details of the abundance distribution 
through emerging bumps in its time evolution. Due to the temporal evolution of the 
thermalization efficiency (see below) and depending on the detailed yields, heavy 
nuclei with p decay half-lives around 14 days can be relevant for the light curve 
behavior at a few weeks, while a few actinides and transuranium elements with 
a-decay half-lives of several tens of days can affect the kilonova emission up to 
hundreds of day after the merger [75, 77]. 

Thermalization efficiency. Thermalization efficiency is defined as the ratio of 
the total energy released by all radioactive processes, nuci, to the energy effectively 
transferred to the ejecta, êth, i.e. fin = étn/énuci- Since both these quantities changes 
with time, also fin is a function of time. Detailed studies on fih can be found, e.g., 
in Refs. [6, 23]. 

First of all, the thermalization efficiency depends on how relevant each decay 
channel is and on how energetic the particles in the final state are. These two points 
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ultimately relate with the physics of the decays and with the actual abundances in the 
ejecta. Beta decay is the dominant decay channel for nuclei A < 200 and it is relevant 
for all Y, conditions, especially for Y, = 0.25. The Q-values of the B decays relevant 
for the kilonova emission are of the order a 1-2 MeV and most of this energy is 
emitted in form of y rays (~ 45%, emitted from the excited daughter nucleus) while 
e~’s takes usually 20%, the rest being lost in neutrinos. The Q-values and the e7 
energies are larger for larger Y.. Alpha decay is the dominant decay channel for 
A = 200. Thus, it affects significantly the heating rate if actinides are produced, i.e. 
if Y, S 0.20. For these decays, typical Q-values are in the range 5-9 MeV. Most of 
this energy is taken by the œ particle while the daughter nucleus de-excitation is 
negligible. Fission is very effective for super-heavy nuclei (A = 250) produced for 
Y, X 0.10. The Q-value of these reactions can be approximated by the kinetic energy 
of the fission fragments and the latter can be estimated as the Coulomb repulsion 
energy of the two fragments at their formation, Qfission ~ Ecou = Z1Z2€7/ (ro (aV oe 
Ay! D; where we have approximated the radius of a nucleus of mass number A as 
roA!/3 with ro ~ 1.8 fm. For typical super-heavy nuclei with (A,Z) ~ (250,100), 
and daughter nuclei with (A ,2,Z1,2) ~ (132,50), Q is ~ 100 MeV. 

Additionally, f} depends on how efficiently the different final state particles ter- 
malize in the plasma. This in turn depends on the physical processes providing ther- 
malization and on the density of the medium, since fin is larger in high density ejecta 
(e.g. in more massive and slower ejecta). High energy charged particles (e~’s, @’s 
and fission fragments) lose their energy in the ambient plasma through Coulomb 
interactions with the free and the atomic e~’s. In the latter case, they can excite or 
even ionize atoms and these are the most efficient thermalization mechanisms. Each 
of these distant interactions transfers a relatively small amount of energy to other 
electrons. Then many interactions are required to thermalize a single projectile par- 
ticle, but due to the low transferred energy the target electrons in the final state ther- 
malize very rapidly. Due to its ZıZ2 dependence, Coulomb interactions with ions 
are relevant only for fission fragments, while strong nuclear interactions with other 
nuclei are negligible. The cumulative nature of the Coulomb processes that affect 
charged particles allows to thermalize efficiently at least a fraction of the available 
energy. Also y-rays lose their energy in the plasma by interacting with electrons 
through photoionization and Compton scattering. Due to the relatively high ioniza- 
tion thresholds of heavy elements (~100 keV), photoionization is the most relevant 
process up to ~ | MeV, while for higher photon energies the Compton scattering 
becomes dominant. 

Any thermalization process is efficient as long the timescale over which it acts 
is smaller than the ejecta expansion timescale. Due to the relatively low opacity 
provided by Compton scattering and photoionization, high energy photons stop to 
thermalize when matter become transparent to them, and for typical kilonova con- 
ditions this happens within the first day after merger (tineff,y ~ 1 day, where tinegt 
is the timescale when thermalization becomes inefficient). For supra-thermal elec- 
trons and @ particles, an efficient thermalization can occur up to several days af- 
ter the merger (tinet ~ 8 day), while fission fragments (again due to the « ZZ» 
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dependence of the Coulomb interaction) thermalize efficiently up to a few weeks 
(finet ~ 16 day). Detailed calculations (e.g. [6]) show that fin(t) decreases from 
0.5 — 0.6 during the first day down to 0.1 around 10 days after merger. 

Atomic opacity. The photon opacity, Ky, quantifies the degree of transparency of 
matter to electromagnetic radiation. In particular, it can be understood as the cross 
section (0) per unit mass of a fluid element to radiation: Ky = no/p, where p is the 
matter density and n the target particle density. In general, Ky depends on the energy 
of the incident photon and on its complex interactions with the electron structure of 
the atom, in all its possible ionization states. Thus information about the composi- 
tion and the ionization degree of each species is crucial to take properly into account 
the most relevant atomic opacities. While the ejecta expand and cool, electrons re- 
combine with atoms to form ions and neutral atoms. First ionization energies vary 
between 3 and 25 eV, while the innermost electrons of heavy elements have ion- 
ization energy of ~ 100keV. Given the expected temperature, Eq.(38), most of the 
electrons have recombined at the time of the kilonova emission. and the abundance 
of free electron is Y; free ~ a few 0.01, decreasing as a function of time due to elec- 
tron recombination. Despite the fact that plasma collisions are not effective enough 
to maintain thermodynamics equilibrium during the ejecta expansion, due to the 
high opacity expected in the ejecta of compact binary mergers, radiation can effec- 
tively drive the ion abundances towards local thermodynamical equilibrium (LTE), 
at least during the first days after merger (e.g. [31] and references therein). 

For the photon energy interval relevant for kilonovae, bound-bound transitions 
are the most important atomic processes, followed by bound-free and free-free 
opacities, usually smaller by several orders of magnitudes. Detailed and exhaus- 
tive experimental values of Ky for heavy elements, for a broad range of photon 
energies and in relevant thermodynamics conditions, are mostly missing, and the 
recourse to sophisticated, but still uncertain, atomic physics calculations is neces- 
sary (e.g. [30, 64]). The energy-independent electron scattering opacity (Thomson’s 
opacity), Ky,rh = 82/3 (ahic/mec)” ne free /p ~ 3.97 x 107? cm? g7! (Ye free /0.01), 
where œ = 1/137 is the electromagnetic coupling constant, sets a lower limit that be- 
comes relevant only for photon wavelengths Ay = 104A. Energy- and composition- 
averaged values of the opacity (sometimes called gray opacity) can be use to roughly 
characterize the global plasma behavior. For example, for the Planck mean opacity 
K,(Ay) is averaged over the Planck distribution function, while for the Rossland 
mean opacity 1/,(Ay) is averaged over the temperature derivative of the Planck 
distribution function. 

When the ejecta contain light r-process elements (A < 140), bound-bound tran- 
sition involve d-shell valence electrons. For Ay = 10°A the opacity strongly de- 
creases with the photon wavelength (Ky(Ay = 10°A) ~ 10°cm? g7! while Ky(y= 
10*A) ~ 10~3cm? g~!) and typical gray opacities are < 1cm? g™!. If lanthanides 
and/or actinides are also present, the opening of the electron f-shell increases enor- 
mously the number of possible transitions and the bound-bound opacity is char- 
acterized by a forest of lines, Doppler-broadened by the large expansion velocity. 
The spectral opacity still decreases with the photon wavelength, but much mildly 
(Ky(10°A) ~ 103cm? g~!, but K,(2.5 x 104A) ~ 107°cm? g~!), and typical gray 
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opacities are in this case > 10 cm? g~!. We stress that even a small amount of lan- 
thanides or actinides (X,4/, +XĘavs < 1074) can already change the global matter 
opacity. 


Modeling kilonovae 


Kilonova modeling is an extremely challenge task. State-of-the-art models require 
the solution of the photon radiative transfer equation in the expanding ejecta, pos- 
sibly in multiple dimensions, e.g. [31, 63, 74]. In addition to physically motivated 
density and temperature profiles, theyse models also require detailed information 
on the nuclear composition, on the ion abundances and on the spectral opacities 
at every time and everywhere in the ejecta. These models are able to provide light 
curves and spectra from a given configuration, at different epochs after merger. If 
the measurement or computation of all relevant transition lines is one of the ma- 
jor theoretical uncertainties, also their treatment in radiative transfer codes and in 
particular their translation in an effective opacity (“effective” because it refers to 
a discretization procedure involving finite wavelength interval, where many single 
lines are present) is not obvious. 

In the following, we will derive fundamental scaling relations, based on a simpli- 
fied analytical model, which are in qualitative agreement with more detailed mod- 
els. In addition to provide the fundamental scales of the problem, they also highlight 
the impact of nuclear physics input and the variety implied by the different ejecta 
conditions expected in compact binary ejecta. To model the ejecta we consider a 
spherically symmetric distribution of total mass mej and average speed (vej), char- 
acterized by a grey opacity Ky and expanding homologously. Matter at the outer edge 
is moving at velocity Ve; max, and at each time f its radial position is Rmax = Vej,max?. 
Any internal shell of mass 6m is expanding at a speed Vej < Vej,max, Constant in 
time and proportional to the radius, such that its radial position evolves accord- 
ing to R = vet. We further define Meny the mass of the envelope above ôm and 
AR = Rmax — R its radial thickness. The expansion timescale of this envelope can 
be computed as fexp ~ AR /vej. Thermal photons produced inside the envelope will 
contribute to the kilonova if fexp ~ taiff, Where taig¢ is the photon diffusion timescale. 
The latter can be determined from random-walk arguments starting from the photon 
optical depth (Ty) and mean free path (ly = 1/ (Kyp)) as taife ~ Tyly/ c. The optical 
depth is defined as the integral of the photon inverse mean path, £y ' along an outgo- 
ing radial path. In words, Ty counts the average number of interactions that a photon 
experiences before being emitted at the photosphere and it can be approximated by 
Ty ~ (p)KAr, where (p) is the average density experienced by the outgoing photon, 
(P) ~ meny/(472R? AR). By equating the expansion and the diffusion timescale we 
can determine the time 7(meny) at which the photons emitted by Meny will contribute 
to the kilonova, and we can estimate the peak time of the kilonova emission fpeax by 
taking meny ~ Mej and vej ~ (Vej): 
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Mej Ky Ky 1/2 Mej 1/2 (vej) —1/2 
4rlvej)c 10 cm? g 0.01 Mo O.1¢ 


(43) 
The energy available to power the kilonova at f is the nulcear energy released by mej 
and thermalized by the plasma: Ly(7) ~ êo (1/1 sec)~* fin(#) mraa(f), with œ = 1.3 
as a typical value. Once again, we can estimate the peak luminosity by taking f = 
tpeak ANd Meny ~ Mej to obtain: 


-a/2 1-a/2 
L ak ~ 2.4 x 10" erg s7! =a “i sue ee 
i 10 cm? g7! 0.01 Mo 


(vej) n éo Sin 
G 3 (srar l =) (2). ae 


The radius of the photosphere at f is Rpn(f) © vejf and then at the luminosity peak: 


1/2 i 1/2 A \ 1/2 

K m (vej) 

Ronneak ~ 1.26 x 1015 i = . 5 
ph,peak 3 = ( 10 cm? =) (oa Mo 0.lc a 


Finally, assuming black body emission and using the Stefan-Boltzmann law, the ef- 


fective photospheric temperature, Ty eff, can be determined as T$ = (Ly / (422, 3p) ) ; 


This expression can be evaluated at the peak: 


3 Ky —(a+2)/8 Whe —a/8 
Tees 2.15 x 10°K ( = SA 
arr (o cm? g! ) (on iz) 


(vej) \ 072/8 ĉo V4 7p yA a 
O.lc¢ 5 x 10!erg s7! g~! 0.5 


and translated in a peak wavelength Ay peak = 1.35 x 10°nm (T; oe 15 x 10°K)~! 

The above estimates have been done assuming Ky ~ 10 cm Le. considedns 
ejecta that contain a significant fraction of lanthanides and sas (meaning that 
the ejecta contains some matter with initial Y, < 0.25). In this case, the luminosity 
peak is expected to occur around one week and at near-IR wavelengths. The lu- 
minosity is more than 100 times larger than the one of a typical nova, but several 
orders of magnitudes lower than a supernova (SN). Since SNe are also powered 
by radioactive decays in expelled matter, the mean reason for such a large differ- 
ence is in the amount of ejecta (Ly « me ay 3, much lower in the case of compact 
binary mergers. If the wind ejecta or the dynamical ejecta have been significantly 
irradiated by neutrinos, the initial Y, could have increased such that the production 
of lanthanides is prevented. In that case, K < 1cm? g7! and the lumonisity peak is 
expected to happen earlier (around 1 day), with a higher luminosity (more than 10° 


times the one of a nova) and at bluer peak frequency (Ay ~ 500 nm). 
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Fig. 11 Modeling of the bolometric light curve of AT2017gfo, the kilonova associated with the 
GW170817 BNS merger event. The slope of the light curved, fitted with four different models 
characterized by different ejecta mass and opacity, is consistent with a power-law decay, t~’ such 
that B ~ 1.0 — 2.0 (@ in the text), compatible with the decay of r-process material in the ejecta of 
compact binary mergers. Figures taken from from [60]. 


GW170817 and its kilonova 


On August the 17th 2017 and in the subsequent weeks, the first unambiguous kilo- 
nova resulting from a compact binary merger was detected, see e.g. [3]. This kilo- 
nova (called AT2017gfo) followed GW170817, the first GW signal compatible with 
the late inspiral of two NSs [2]. The total mass of the system was 2.74 Mowhile the 
mass ratio, assuming slowly spinning NSs, was measured to be between 0.7 and 1.0. 
Light curves in different photometric bands of this unprecedented UV/visible/IR 
emission showed an early peak (around 1 day after merger) in the visible frequen- 
cies, followed by a later peak (around 5-7 days after merger) in the near-IR, see 
[68] and references therein. The bolometric luminosity of the event 1.5 days after 
merger was ~ 3.2 x 10*erg s~!, while it decayed approximately following a £713 
power-law during the first week so that around 7 days Ly ~ 6 x 10*°erg s~!. The 
spectrum at 1.5 days was very close to a black-body of 5 x 10°K, while around 7 
days it was broadly compatible with a 2.2 x 10°K black-body spectrum with a forest 
of absorption features, e.g. [49, 60, 65]. The identification of elements in the spec- 
trum is very challenging due to the high density of lines and to their broadening due 
to the high expansion velocities. An analysis of the spectra recorded during the first 
days has revealed features compatible with the presence of Sr, an element of the first 
r-process peak [72]. 

Theoretical modelling of AT2017gfo requires the presence of more than one 
component of the ejecta. The different components are characterized by different 
masses, velocities and opacities, and possibly a non-trivial dependence from the ge- 
ometry of the ejection. This is indeed necessary to explain the observed blue and red 
peaks [68, 48]. For example, one can consider the presence of two distinct compo- 
nents to explain the color evolution of AT2917gfo. Using Eq. (43)-(46), assuming 
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Fig. 12 Elemental abundance comparison between BNS merger models and metal poor star obser- 
vations. The violet curve represents nucleosynthesis yields for tidally dominated dynamical ejecta, 
while the other curves yields from neutrino-driven wind ejecta for different massive NS lifetimes. 
Theoretical abundances are compared with two classes of metal poor stars. In the right panel dy- 
namical ejecta have been diluted by a factor of 50 with respect to the neutrino-driven wind one. 
Figures taken from [40]. 


a œ ~ 1.3 and a thermalization efficiency of 0.8 and 0.4 for the peaks at 1.5 and 7 
days, respectively, the above peak properties imply that the blue peak was character- 
ized by Mpjue © 0.019 Mo, Ky blue © 1 cm? g7! and vej blue © 0.20c, while for the red 
peak mrea © 0.058 Mo, Ky,rea © 4.2 cm? g7! and Vejred % 0.09c. The estimated total 
amount of ejecta is thus of the order of several percents of a solar mass. Results from 
one of these models are presented in Figure 11. According to more detailed models, 
e.g [31, 63, 74], the amount of ejecta in this event was ~ 0.02 — 0.05 Mọ. Nuclear 
physics input (e.g. the nuclear mass model) can introduce an additional uncertainty 
factor of a few, possibly reducing the total ejecta mass. The inferred opacities sug- 
gest a negligible amount of lanthanides in the blue component and a lanthanide mass 
fraction between 107? and 107? in the red one. 

The emerging picture is certainly very compatible with results of compact binary 
merger simulations and nucleosynthesis calculations. In particular, the decline rate 
of the bolometric light curve and the inferred opacities are compatible with what 
expected from the collective decay of freshly synthetized r-process elements (see 
for example Figure 11). 


Compact binary mergers and the chemical evolution 


Due to the production and ejection of heavy nuclei, compact binary mergers are a 
possible key player in shaping the evolution of chemical abundances in the Universe. 
A first relevant question is whether they are primary astrophysical sites for the pro- 
duction of r-process elements. Two main ingredients are necessary to answer this 
question: detailed abundance predictions and reliable merger rates. The former are 
obtained by combining the amount of ejecta predicted by merger simulations with 
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Fig. 13 Cumulative mass fraction obtained by the solar residual r-process abundances from [51] 
both as a function of the mass numbers (left panel) and the atomic numbers (right panel). 


the distributions of the nucleosynthesis yields corresponding to the properties of the 
ejecta. The latter are related to the binary formation channels and evolution. In par- 
ticular, they depend on the probability that a compact binary forms and on the time 
that it takes for it to reach coalescence due to GW emission. To provide a simple an- 
swer, we compare the amount of r-process elements in the Milky-Way (MW, M,.mw) 
to the whole amount of r-process material ejected by compact binary mergers dur- 
ing the galactic evolution (M;.cgm). Their ratio, fcgm = M,,CBM /M;mw, quantifies 
the relevance of compact binary mergers in accounting for the observed r-process 
elements. The enumerator can be evaluated as the product of the average merger 
rate in the Galaxy, Rcgm, times the average mass of r-process elements expelled by 
a single event, m,.cpm, times the age of the Galaxy, tmw + 13.5 Gyrs. The mass per 
event can be estimated from the lower bound obtained from GW170817 ejecta, i.e. 
m,cpm ~ 0.02M.. The merger rate is still very uncertain, but it can be measured in 
different ways. From the GW events detected so far in the first two observing runs, 
the Ligo-Virgo collaborations have provided a rate of 110 —3840Gpc~? yrs! for 
BNS mergers, and a robust upper bound of 660 Gpc~? yrs~! for BHNS mergers 
[5]. Assuming that a significant fraction of mergers produce also a short gamma- 
ray burst (SGRB), the merger rate can be inferred also from the SGRB rate. The 
latter is bracket by Rsgrg = 0.6 — 6 Gpe~*yrs! [71, 18], and since GRBs are colli- 
mated emissions with a beaming correction factor fbeam ~ 100 (fbeam = 1 — cos Giet 
and typical opening angles are jet ~ 5 — 10°, [15]), the resulting lower bound on 
Rcpon rate is ~ 80 —600 Gpe~3yrs~!, More theoretical bounds can be obtained from 
population synthesis studies and from extrapolations of the observed populations 
of galactic BNS system (see for example [1] and references therein). These rates 
broadly agree with the measured rates, even if with larger uncertainties. Assuming 
a galaxy density of ~ 0.01 galaxy Mpc~?, a conservative merger rate inferred from 
theoretical models and observations is Regm % 30 Myr~!galaxy~!. For the total 
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mass of r-process elements in the MW, we consider the Sun as representative of the 
averaged enrichment of the Galaxy, a common assumption in Galactic chemical evo- 
lution models. We also assume that all this mass is produced by the one single type 
of events and with a relative constant yield and rate for the entire evolution of our 
Galaxy. These hypothesis are also grounded in the observations of the abundances 
of r-process elements in metal poor stars. Indeed the spectroscopic analysis of the 
light emitted by a star reveals the chemical composition of the ISM from which 
the star has formed. Low mass stars can live longer than the Galaxy, thus some of 
them are among the oldest objects in the MW. These old stars have an extremely 
low metallicity since they formed very early in the Galactic history when only a few 
SNe had exploded and polluted the ISM. They can be considered fossils of the early 
chemical enrichment of our Galaxy and the chemical abundances measured in their 
spectra can be studied to infer the characteristics of the first stellar nucleosynthesis 
events. Noteworthy, in the spectra of these old stars the lines of neutron capture ele- 
ments are identified and their abundances measured. Among these fossil stars, there 
is a group that appears particularly rich in r-process elements (so called r-process 
rich stars, see e.g. [61]). For this group, the derived abundance pattern reproduces 
very closely the solar residual r-process pattern, at least beyond the second peak, 
while the first peak presents a significant dispersion. Since the low metallicity im- 
plies that these stars have been polluted by very few (if not one single) r-process 
sources, the common patterns observed in the solar r-process residual, in this group 
of metal poor stars and in compact binary merger nucleosynthesis calculations re- 
veal the presence of robust features that characterize the r-process nucleosynthesis 
in very general terms. Nevertheless, the emerging picture is more complicated than 
that: other metal poor stars present a substantial amount of first peak elements and a 
lower enrichment in the heavy ones [22]. This is still compatible with the large vari- 
ability expected in the the yields of compact mergers (see Figure 12), but it leaves 
also space for other production sites. To account for these uncertainties, starting 
from the solar r-process abundances, we compute the cumulative mass fraction of 
the r-process materials above a certain mass number Amin, X> (A > Amin), see Fig- 
ure 13. To estimate the typical total mass fraction of r-process elements in the MW 
we assume X;— proc © X> (A > Amin) and we consider Amin = 68,89, 124 such that 
X;—proc(A > Amin) = 40,9, 6 x 1078, respectively. Since the mass in stars and ISM 
of the MW is Myw 7 6 x 10!°M 5, Scop can be finally estimated as: 


RcsM Mr, CBM tMw X;—proc 
~ 1.35 47 
fesm (ses) Ge (= | ( 10-7 ) Eo &) 


and fcgm ~ 0.34, 1.50,2.25 for Amin = 68,89, 124, respectively. 

These estimates are rather crude in many ways. For example, they assume that 
the present day merger rate is representative of the average merger rate, while we 
know that it has significantly changed during the Galaxy history. Chemical evolu- 
tion models can be used to better consider the merger rate and set more stringent 
constrain on the impact of the enrichment of compact mergers to the total balance 
of r-process in the Galaxy, see [11]. Moreover, it is not obvious that the reference 
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Fig. 14 Left: Scatter plot of the relative abundance ratio of iron over hydrogen versus europium 
over iron, for a large sample of metal poor stars. The quantity [A/B] is relative with respect to 
the solar one, i.e. it is defined as [A/B] = log; (Ya/Yp) — logio (YA /YB)o. Data have been taken 
from [56]. Right: Summary of the observational evidences for the production of r-process elements 
expressed on the r-process mass (m,) versus galactic event rate (Rmw) plane. For concreteness, r- 
process abundances are assumed to be equal to the solar residual for Amin = 90 while the BNS rate 
is inferred from galactic BNS observations [50]. (Courtesy of K. Hotokezaka; see also [24]). 


amount of ejecta is representative for the merger population. However, it is worth 
noticing that an IR excess, observed in the afterglow light curves of a few SGRB, is 
broadly compatible with the emission expected from a kilonova and it usually points 
to the presence of large ejecta masses (107? — 107! Mo), e.g. [66]. In addition, as- 
suming that a significant fraction of the accretion disk is expelled in disk winds, our 
reference value is also well compatible with the results reported in Figure 3. Some 
astrophysical processes are also neglected for simplicity. For example, when com- 
pact objects emerge from CCSNe, they receive a kick at birth. This kick, in addition 
of being a threat for the survival of the binary, can put the binary in a wide orbit 
inside the galactic potential, such that the merger could happen with a significant 
offset from the stellar and gas distributions inside the Galaxy (see [24] for a more 
detailed discussion). Despite the large uncertainties, fcgm ~ | and this testifies that 
compact binary mergers are primary astrophysical site where r-process elements 
are produced. This is especially true for elements above the second r-process peak, 
while the explanation of the first peak could be a clear confirmation that different 
mergers can produce different yields or could require additional production sites. 

A related question is whether other astrophysical sites are able to synthesise 
heavy elements through the r-process. Historically, proto-neutron star winds emerg- 
ing after a successful CCSN explosion have long been thought to be a possible 
r-process nucleosynthesis site [52]. The argument applied above to relate the rates 
and the ejecta from compact binary mergers to the observed amount of r-process 
elements in the MW could be easily adapted to regular CCSNe and it is not able 
to distinguish alone between these two scenarios. However CCSNe are signifi- 
cantly more frequent than compact binary mergers and their rate is well constrained, 
Rccsyn % 2.8 x 104Myr"! (this rate can be related to Rcgm by considering that ap- 


r-process nucleosynthesis from compact binary mergers 47 


proximately half of the massive stars are in binaries, two CCSNe are required to 
form a compact binary and only ~ 1% of the stellar binaries survives the two CCSN 
explosions). From Eq.(47) it is evident that CCSNe could explain the bulk of the 
galactic r-process nucleosynthesis if every SN would eject ~ 1075 Mo of r-process 
elements. Thus, the competition here is between rare events that expel large amounts 
of r-process material (e.g. compact binary mergers) or frequent events with much 
smaller amounts (e.g. regular CCSN). Over the past years, several observational ev- 
idences (in addition to GW170817) have accumulated pointing to the fact that the 
r-process elements come from rare events that produce a significantly large amount 
of r-process elements. In the following we will briefly review them: 

Eu abundance in galactic metal poor star. The analysis of the abundance of r- 
process elements, and especially of Eu, as a function of metallicity (and in par- 
ticular of iron and of œ elements, which are good tracers of the early evolution 
of the ISM composition due to the explosions of the very first massive stars) has 
revealed that the average abundance of Eu correlates with the one of œ elements 
produced in CCSNe through the entire metallicity evolution (it correlates also with 
the one of iron, but only before type Ia SN start to explode, i.e. for low metallicity). 
However, differently from the œ elements case, the distribution of the single obser- 
vations changes considerably: while at present metallicity the cumulative effect of 
several nucleosynthesis episodes and the efficient mixing of gas inside the Galaxy 
has homogenized the ISM composition and reduced the observed spread, at early 
times, when the amount of iron was 10% — 10? times smaller than the present one, 
the ratio of the Eu over iron abundances shows a two orders of magnitude scatter, 
ranging from metal poor stars where Eu is underrepresented to cases where this 
ratio is more than 10 times the one observed in the present solar system (see the 
left panel of Figure 14 and [61]). This large scatter suggests that the possible r- 
process elements pollution comes from a single rare event that ejects large amount 
of r-process material in ISM clouds that have seen only a few SNe. Considering 
that ISM mixing over the entire Galaxy requires a timescale much larger than the 
star formation timescale, this naturally introduces the inhomogeneous character to 
the early chemical evolution necessary to explain the observations. Sophisticated 
Galactic chemical evolution models have been specifically developed to mimic this 
observational spread (see [1 1] and references therein). It shall be underlined that this 
is a peculiar characteristics of the r-process elements, not present in other chemical 
elements, such as @—elements produced together with iron by the more frequent 
SNe. 

r-process abundance in ultra-faint dwarf galaxies. In addition to galactic ob- 
servation, Eu can be observed also in classical dwarf and ultra-faint dwarf (UFD) 
galaxies, satellites of the MW and formed by very old stars (interestingly, these ob- 
servations are possibly related with the ones of galactic metal poor stars, because the 
latter are often located in the so called galactic halo, that is thought to be formed by 
accretion episodes of dwarf galaxies on the MW). In particular, in the case of UFD 
galaxies, while Fe is observed in all cases, Eu has been detected so far only in a cou- 
ple of cases (Reticulum II and Tucana III), corresponding to ~30% of the available 
sample [28, 20]. Since UFD galaxies are formed by 10° — 10° stars, such a low de- 
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tection rate points again to a single, rare event (no more than one per 104 — 10° stars). 
Moreover, it is interesting to compare the ratio between the amount of iron and of 
europium observed in Reticulum II (Mpe Ret n ~ 0.7 Moand Meu ret 1 © 1075Mọo), 
with the theoretical value predicted by assuming that Fe is produced by CCSNe 
(typically, mfe,ccsn © 0.07Mo) and Eu by merger events (for which we conserva- 
tively assume again m;cpm ~ 0.02Mo) with an abundance equal to the solar one 
(XEu © 4.1 x 10710). Since (Mre /Mtu)iheory ~ MFe,ccsNRccsn/MEu,cBMRcBM and 
MẸu,CBM ~ Mr.CBM XEu/Xr-proc» We finally obtain: 


(Mre /MEu)ineory 
(Mre /MEu)Ret 1 


= 5.69 poate as (Zenit) 


7x 104 4.7 x 102 


MFe,CCSN /Mr,CBM XEu aE (48) 
7/2 4.1 x 107210 10-7 } ` 


Thus, despite the crude approximations, the large uncertainties and the possible pe- 
culiar character of Reticulum II, the relative amount of Fe over Eu is broadly com- 
patible with what expected from compact binary enrichment. 
Radioactive elements in Earth sediments and in meteorites. The analysis of the 
abundances of long-lived r-process isotopes (¢)/2 > 10’yr) and of their daughter 
nuclei in meteorites (formed when the solar system was formed) and in sediments 
(e.g. deep see floor) reveals information on the local isotopes production at specific 
times (see [24] and references therein). Due to the paucity of the r-process nuclides 
129] and 747Cm, whose lifetimes are < 20 Myr, and to the relative large abundance 
ratio (0.008) between 7**Pt (ti j2~ 80 Myr) and its stable daughter nucleus 238 
in meteorites, it is possible to set the delay time between the solar system formation 
and its last r-process pollution to ~ 100 — 120 Myr. In addition, the analysis of the 
deposition rates in deep see sediments of radioactive nuclides like Fe (very short- 
lived nuclide usually produced in SNe) and the r-only 7“*Pt over the past 25 Myr 
shows a large fluctuation, with the present rate ~ 10~* times smaller than the rate 25 
Myr ago. This indicates that the enrichment of r-process elements is uncorrelated 
with the CCSNe enrichment and the former is due to rarer events that produce large 
yields [25]. 

A summary of the above observational constraints is provided in the right panel 
of Figure (14). All of them work against regular CCSNe and in favor of compact 
binary merger as major site for the r-process production. A final interesting question 
is whether compact binary mergers are the only sources of r-process elements. The 
large abundances of the first peak elements in the solar r-process pattern and in some 
metal poor stars, with respect to the elements beyond it, could be an indication of 
the occurrence of mergers in which the production of heavy r-process elements is 
disfavored (for example, due to a strong neutrino irradiation). However, it leaves 
also space for other sources able to provide weak r-process nucleosynthesis (e.g. 
electron capture SN, [69]). Moreover, it is non-trivial for compact binary mergers 
to explain the abundances at very low metallicity. Indeed, a compact binary merger 
requires the successful explosion of at least one CCSN to form it (two for a BNS 
system) and, in addition to the stellar evolution timescale, it is necessary to wait for 
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the GW-driven inspiral timescale. The latter can be estimated from the semi-major 
axis, a, the eccentricity, e, the total and reduced mass, M and u, of the binary at 
formation as: 


a 4 M = u -1 2\7/2 
taw ~ 0.66 Gyr 1 . (4 
eae o. (53 i) (z) one) i 


This is a possibly rather long timescale to be reconciled with fast mergers happening 
at extremely low metallicity. However, the strong dependence of tgw on a and e does 
not exclude them, but it requires a population of tight, possibly eccentric, compact 
binaries. Possible alternative sites are represented by rare classes of SNe character- 
ized, for example, by very intense magnetic fields and fast rotating cores (e.g. [73]). 
It is important to stress that the relation between time and metallicity is not obvious: 
yields expelled by stars, SNe and compact binary mergers need to mix with the ISM 
before entering in the composition of the next stellar generations. Moreover, due 
to kicks and fast traveling ejecta, the places where SNs explode, compact binaries 
merge and new stars form could be distinct. All these effects are amplified in low 
metallicity conditions, where only a few enrichment episodes have happened. 


Summary and outlook 


Compact binary mergers involving at least one neutron star represent ideal environ- 
ments where the production of heavy elements through the r-process nucleosynthe- 
sis takes place. The prediction of the precise composition of the ejecta and a clear 
understanding of its origin depend both on the modeling of the astrophysical sites 
and on detailed nuclear physics knowledge. A lot of progress has been achieved 
in the past few years and the predicted abundances are able to explain many inde- 
pendent observations. Compact binary mergers are also possibly very relevant in 
explaining the nuclear evolution of the Universe, in terms of the abundances that we 
observe at different epochs and in different astrophysical environments. 

Still, many open questions remain. On the one hand, more realistic and sophisti- 
cated compact binary merger models are required to predict the properties and the 
amount of ejecta with better accuracy and properly taking into account all the rele- 
vant physics. On the other hand, a more robust knowledge of the properties of exotic 
neutron-rich nuclei is key to reduce present nuclear uncertainties. In this respect, 
existing and upcoming world-wide nuclear facilities (including FAIR, FRIB, HIAF, 
RAON, RIKEN and SPIRAL) will finally produce some of the neutron-rich nuclei 
relevant for the r-process and measure their properties. Multimessenger observa- 
tions, as well as the study of the composition of matter in different astrophysical 
and terrestrial contexts, will also sharpen our understanding and test our models, 
helping reducing our ignorance and forcing us to look at the problem from many 
viewing angles. Despite a conclusive answer about the presence and the role of 
other possible r-process nucleosynthesis sites still needs deeper investigations and 
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clear evidences, it is nowadays certain that compact binary mergers are one of the 
major sources of r-process elements in the Universe. 
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